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We report here results from a Chandra ACIS observation of the stellar populations in and 
around the Messier 17 Hn region. The field reveals 886 sources with observed X-ray luminosities 



' (uncorrected for absorption) between ~ 29.3 < \ogL x < 32.8 ergs s 1 , 771 of which have stellar 

^"*^ ■ counterparts in infrared images. In addition to comprehensive tables of X-ray source properties, 



several results are presented: 

1. TheX -ray Luminosity Function is calibrated to that of the Orion Nebula Cluster population 
to infer a total population of roughly 8000-10,000 stars in M17, one-third lying in the central 
NGC 6618 cluster. 

■ 2. About 40% of the ACIS sources are heavily obscured with Ay > 10 mag. Some are concen- 

ly-j , trated around well-studied star- forming regions — IRS 5/UC1, the Kleinmann- Wright Ob- 

ject, and M17-North — but most are distributed across the field. As previously shown, star 
formation appears to be widely distributed in the molecular clouds. X-ray emission is de- 
f~>) \ tected from 64 of the hundreds of Class I protostar candidates that can be identified by 

near- and mid-infrared colors. These constitute the most likely protostar candidates known 
in M17. 

q • 3. The spatial distribution of X-ray stars is complex: in addition to the central NGC 6618 

cluster and well-known embedded groups, we find a new embedded cluster (designated 
M17-X), a 2 pc-long arc of young stars along the southwest edge of the M17 Hn region, 
and 0.1 pc substructure within various populations. These structures may indicate that the 
populations are dynamically young. 

4. All (14/14) of the known O stars but only about half (19/34) of the known B0-B3 stars in 
the M17 field are detected. These stars exhibit the long-reported correlation between X-ray 
and bolometric luminosities of L x ~ 10~ 7 L}, o i. While many O and early B stars show the 
soft X-ray emission expected from microshocks in their winds or moderately hard emission 
that could be caused by magnetically channeled wind shocks, six of these stars exhibit very 
hard thermal plasma components (kT > 4 keV) that may be due to colliding wind binaries. 
More than 100 candidate new OB stars are found, including 28 X-ray detected intermediate- 
and high-mass protostar candidates with infrared excesses. 
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5. Only a small fraction (perhaps 10%) of X-ray selected high- and intermediate-mass stars 
exhibit K-band emitting protoplanetary disks, providing further evidence that inner disks 
evolve very rapidly around more massive stars. 

Subject headings: open clusters and associations: individual (M17) — X-rays: individual (M17) 
— stars: early-type — stars: pre-main-sequence — X-Rays: stars 



Introduction 



In the last 30 years, most observational studies of the stellar denizens of massive star-forming regions 
(MSFRs) were based on broadband visual photometry, supplemented by occasional visual spectroscopy of 
the most massive cluster members, and near-infrared (NIR) photometry of small regions (< 10') situated 
on the cluster cores. In a few visually obscured regions this was augmented by NIR spectroscopy of the 
brightest OB cluster members. The surrounding Hn regions and molecular clouds were often well-studied in 
radio/millimeter wavelengths, including arcminute-scale maps of molecular line emission and arcsecond-scale 
maps of continuum structures. More recently, heated dust structures have been mapp ed in the mid-infrared 
(MIR) with space-based instruments on the Midcourse Spac e Experiment (MS X) (e.g. iKraemer et al.l 120031 ) 
and the Spitzer Space Telescope (e.g. IChurchwell et al.ll2004t lAllen et al.ll2005h . 



In the high-energy regime, lower-mass young stars in nearby star-forming regions were studied ex- 
tensively with X-ray missions during the 1980-90s ([Feigelson fc Montmerldll999| ). but studies of the more 
distant MSFRs were hampe red by the limited sp atial resolution of these telescopes. The launch of the 
Chandra X-ray Observatory (jWeisskopf et al.l 120021 ) in 1999 makes X-ray studies of MSFRs much more fea- 
sible, due to the sub-arcsecond spatial resolution of its mirrors, their high reflectivity at shorter wavelengths 
that penetrate interstellar absorption, and the excellent noise characteristics of its primary imaging camera, 
the Advanced CCD Imaging Spectrometer (ACIS). High-resolution X-ray images of MSFRs can reveal a 
complex menagerie of hundreds of stellar sources (protostars, pre-main sequence stars, and OB/ Wolf-Rayet 
stars) and diff use structures such as wind-swept bubbles, superbubbles, chimneys, and supernova remnants 
(see review bv lFeigelson et al.l 12006). X-ray selected samples are surprisingly advantageous for uncovering 
the stellar populations of young clusters. By studying a number of these fields with varying ages and stellar 
content with Chandra, we hope to characterize these stellar populations and to learn more about the origins 
of any diffuse X- ray emission that m ay be presen t. Some of the m ost recent C handra MSFR studies in- 
clude NGC 6193 Jskinner et al.ll2005l). N GC 6334 fezoe et alJbqod). 3 Doradus JTownslev et al.ll2006allbh, 
Cepheus B/OB3b dGetm an et alJl2006ah. RO W 38 (|Wolk et alJl200d). Westerlun d 1 dSkinner et alJ l200fi 



Muno et al.ll2006bl) . W 49A (jTsuiimoto et all 120061) . and NGC 6357 (|Wang et all 120061) . These studies all 



show rich stellar populations, many uniquely revealed by the high-resolution Chandra observations. 

These X-ray observations are complementary to high-resolution IR imaging data from 2MASS, Spitzer, 
and modern ground-based facilities. We have shown that Chandra observations of Galactic MSFRs are not 
strongly contaminated by X-ray sources that are not cluster members; in a typical Chandra observation of 
a cluster around 1 — 2 kpc distant, simulations and data show that only a handful of foreground stars and 
background active gal actic nuclei (AGN) are present (~5% of the X-r ay sources), while hundreds of cluster 
members are revealed (jGetman et al.ll2005aL l2006at IWang et al.ll2006l ). These X-ray detections can be used 
to identify cluster members for further IR study without the usual selection by IR excess. This provides a 
disk-unbiased sample for further IR study, at least for lower-mass T Tauri (pre-main sequence) stars. X-ray 
observations also are not hampered by the confusion due to bright dust and gas emission that plagues IR 
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and visual studies of Hn regions. X-rays are often the only way to detect fainter cluster members along 
ionization fronts or in the vicinity of bright OB stars. 

Here we present the first high-spatial-resolution X-ray point source study of M17 (<^> W 38 S 45 
RCW 160 <^ the Omega Nebula), a famous and spectacular nearby MSFR (Figure [TJ. M17 is a blister Hn 
region at the northeast edge of o ne of the largest giant molecular clouds in the Galaxy, with an extent of 4° 
(~110 pc. lElmegreen et al.lll979| ). The local geometry rese mbles the Orion Nebula Hn region except that it 
is viewed edge-on rather than face-on ( Meixner et al. 1992). Molecular material is concentrated in a clumpy 
structure known as M17-SW. Active star formation is concentrated in three molecular cores (Figure [T]a) : 
around the ultracompact Hn region M17-UC1 and its assoc iated masers, the double B star known as the 
Kleinmann- Wright Object (KWO) and its associated cluster (|Chini et al.ll2004al ) . and the protostars in M17- 
North. Historically, their proximity to the Hn region was attributed to sequential t riggering of st a r form ation 
(|Elmegreen fc Ladalll977l ). although this link has been challenged for M17-SW bv lWilson et al.l (|2003l ). who 
note that the KWO is too far from the OB cluster and too evolved for its formation to have been caused by 
the cluster. The clumpy molecular structure of the cloud (Figure UJ)) and infrared protostars indicate that 
star formation is still in pr ogress all around the Hn region ( Jiang et al. 20021 Wilson et al. 2003 : Chini et al 
200,4 iReid fc WilsoiJEoolh . 



NGC 6618, the central cluster that illuminates M17, has 100 stars earlier than B9 (jLada et al.lll991 ): 
for comparison, Orion has a dozen. Several obscured 04-05 stars fo rm a central 1' ring and are principally 
responsible for ionizing the nebula ( Hanson. Howarth. fc Contil 1997 ). NGC 6618 is estimated to be no more 
than 1 Myr old; this extreme youth and the absence of evolved stars in the region give good circumstantial 
evidence that we are seeing this cluster before its first supernova, making M17 an ideal place to study the 
interactions of massive stars with their environment and stellar disks before they are affected by nearby 
supernovae. We give a more extensive description of past observations of M17 and NGC 6618 in our earlier 
analysis of th is Chandra observa tion, which focused on the diffuse X-ray emission produced by the cluster 
OB winds (jTownslev et al.ll2003l . henceforth TFM03). 

Distance estimates for M 17 range from 1.3 kpc (Hanson. Howarth. fc Contlll997 ) to 2.2 kpc ( Chini. Elsaesser. fc Neckel 

19801 ) . A thorough study by iNielbock et al.l (|2001l ) gives a dist ance of 1.6 ± 0.3 k pc; we adopt this value 



here. This is consistent with TFM03 and with the NIR study of [Jiang, et all I2002L whose data are used in 
this paper. 

The present Chandra study is based on a 40 ks pointing with the ACIS camera in March 2002; this is 
the same dataset analyzed by TFM03. This image locates over 800 stellar X-ray sources in the NGC 6618 
cluster and cloud population with sub-arcsecond accuracy (Figure [2|) . Chandra detects all of the known 
O stars and about half the cataloged early B stars in NGC 6618, and adds dozens of likely new massive, 
intermediate-mass, and protostellar members as well as hundreds of T Tauri stars. Sparse clusters around 
M17-UC1, the KWO, and M17-North are found. Soft diffuse X-ray emission arising from thermalized O-star 
winds pervades the Hn region, as discussed by TFM03. 

This paper begins with a description of the X-ray sources and their stellar counterparts ($2]). We 
quantify the stellar population and its spatial structure in Two important subpopulations, the embedded 
stars and the massive OB stars, are discussed in detail in the next two sections (SKS|). We summarize our 
findings in fj6] 
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Fig. 1— (a) 2MASS K-band mosaic of M17 covering the 17' x 17' Chandra field. The 115 small diamonds 
mark ACIS source s without identifie d counterparts f £)2.4.ip . Two polygons, reproduced with their labels 
from Figure 2c in Ijiang et al.l (|2002l ) , outline the North Bar and South Bar in the Hn region. Two cir- 
cles, reproduced with their labels from the same figure, outline the central NGC 6618 stellar cluster core 
and the embedded M 17-SW region. A third circle outlines the embedded M17-North region described by 
Henning et al. ( 1998h . Three small boxes indicate the expanded views in Figures fTOl and [TT1 (b) C 18 
contours (-5, -2, 2, 4, 10, 20, 30, 40, 50, and 60 K km s" 1 ) show n on an image from the 2MASS survey, 
taken from the right-hand panel in Figure 1 of IWilson et alj 1|2003[) and reproduced here with permission of 
the authors. This field of view is marked by a large rectangle in panel a. 
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Fig. 2. — A Chandra ACIS-I binned image of M17 with 886 source extraction regions designed to match 
the Chandra Point Spread Function (PSF) shown in blue. The PSF is a complex function of position and 
photon energy; for example the 90% encircled energy radius at 1.5 keV varies from ^0.9" on axis to ~6" at 
the field edge (see §4 of The Chandra Proposers' Observatory Guide at http://asc.harvard.edu/proposer/). 
Contours (green) of bright dust emission from Spitzer 3.6^m data outline the North Bar and South Bar in 
the Hn region. The circles and squares are reproduced from Figure [TJ 
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2. X-ray Observation, Data Analysis, and Initial Results 



M17 was observed with the Imaging Array of the Chandra/ ACIS camera (ACIS-I, 17' x 17' field of view) 
for roughly 40 ks in March 2002. The aim point was centered on the source with the earliest spectral type in 
NGC 6618's centr al ring of O stars, the 04+04 visual binary f £|5.2|) called "Kleinmann's Anonymous Star" 
|Kleinmannlll973h . 



Our data reduction methods are described in detail in Appendix B of TFM03; the same reduced dataset 
used in that study was also used for this one, so we refer the reader to that description for details of the dataset 
and our data analysis. Briefly, the data were corrected for charge transfer inefficiency (CTI) problems in 
the CCDs using the Penn State CTI corrector (jTownslev et al.ll2002h : the 0.5" event position randomization 
added in the Chandra X-ray Center's standard data processing was removed; a sub-pixel positioning scheme 
(jTsunemi et al.l 1200 111 was applied to improve source positions. We then created twelve different images of 
the ACIS field: soft (0.5-2 keV), hard (2-8 keV), and full (0.5-8 keV) X-ray wavebands with four different 
pixel binning scales (4x, 2x, lx, and 0.5x the nominal 0.5" ACIS sky pixel) to cover the full ACIS field, 
the ACIS Imaging Array, the central 8', and the central 4', respectively. 



2.1. Source Detection and Extraction 

A list of candidate point sources was the n obtained in each o f these band-limited images using the 



wavdetect wavelet source detection algorithm ([Freeman et al.l 120021) . The algorithm was purposefully run 
with a low threshold (P = 10~ 5 ) that is sensitive to very faint sources but also identifies some spurious noise 
features. The 10 -5 wavdetect threshold is used because experience has shown that more stringent thresholds 
(e.g. the more typical 10 -6 ) result in many missed sources in regions with complex backgrounds, such as 
M17's diffuse X-ray emission. These twelve source lists were merged (keeping the source position from the 
highest-resolution image) to generate the master list of candidate sources for the full observation. The field 
was visually examined and a few possible faint sources missed by the algorithm were added, resulting in 933 
candidate source positions. 

Photons were extracted using the ACIS Ex tracl[j] (AE) sof t ware p ackage ( Broos et al. 2002h . The pro- 



cedures used in AE are described in TFM03 and lGetman et alJd2005bh . The source extraction regions were 
usually defined to be ~90% contours of the local Point Spread Function (PSF) at 1 keV; crowded sources 
were assigned smaller extraction regions following visual review. AE applies an energy-dependent correc- 
tion to each source's calibration, via the Ancillary Response Function (ARF) file, to account for its finite 
extraction region. 

The background event data set from which local background spectra were constructed was obtained by 
masking regions around all the sources in the catalog via a process that involved two passes through AE. 
First, conservative circular masks covering > 99% of each source PSF were constructed, local backgrounds 
were extracted, and source fluxes were estimated. Then, each source mask was redefined using the generally 
less restrictive criteria that: (a) the mask shall include the source extraction region (polygon) itself; and 
(b) the mask shall also include the region where the expected surface brightness from the point source is 
larger than one half the observed smoothed local background. A circular background region was then defined 



1 The ACIS Extract software package and User's Guide have been available online at 
http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html since February 2003. 
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independently for each source such that the region enclosed ^100 events in the source-free masked data set. 



For each candidate source AE computed the quantity P B , the Poisson probability associated with the 
"null hypothesis" that no source exists with respect to the local background level. We chose two thresholds 
on P B to define the final source list: P B < 0.003 for 845 "primary" sources and 0.003 < P B < 0.01 for 41 
"tentative" sources. Counterparts were identified for 31 of these tentative sources f £|2.4[) , adding confidence 
that these sources are real. The remaining 47 candidate sources had P B > 0.01 and were rejected as likely 
background fluctuations. 

These 886 sources are shown superposed on the ACIST field in Figure [2] AE and the Chandra software 
package CIAO (version 3.2) were used to compute source and local background spectra, calibration filefl 
background-corrected median energies, light curves, variability estimates, broad-band photometry, and other 
quantitative information for each source. Basic properties of the primary and tentative sources are presented 
in Tables [T] and O Sources are listed by increasing right ascension and can be identified by their sequence 
number (e.g. #358) or their IAU desi gnation (e.g. CXQU J1 8202 5. 70-161649.9) Please see the table notes 
for definitions of the columns and see iGetman et al.l (|2005br ) and iTownslev et al.l (|2006br ) for further details 
on the derivations of the properties. The positions reported are estimated by AE; the Chandra coordinate 
system was aligned to the Hipparcos reference frame using 51 matches between strong on-axis ACIS sources 
and reliable 2MASS stars fi j2.4|) . The formal uncertainty in the remaining systematic reference frame offset 
between Chandra and 2MASS is ~0.04" in each axis. 



2.2. Source Variability 



Thirty-nine of the M17 ACIS sources show significant X-ray variability (Pks < 0.005 in column 15 
of Table [l}. We display some of these variable light curves in Figure[3] Many show t he high-amplitude , 
fast-rise-slow-decay morphology characteristic of solar- type magnetic reconnection flares (jReale et al.ll2001 ). 
Some show secondary flares superposed on the decays of primary flares. Others show a poorly understood 
slow rise over several hou rs. This wide range of flare morphologies is characteristic of X-ray luminous pre- 
mam sequence stars (e.g. Ilmanishi et al.l[200lt iGrosso et alJl2004t iFavata et alJl2005t Iwolk et ZlEjoj ), A 
few sources show more random variability not characteristic of T Tauri flares. One bright variable star 
(#396) is not shown here, but is discussed in Section [5731 light curves for 24 variable but weak sources are 
not shown. 



2.3. X-ray Fluxes and Luminosities 

Characterization of the spectral properties and luminosities of the sources in this data set is a challenging 
endeavor because the quantity of data available for most sources is so low: 95% of sources have fewer than 
100 net counts. Such data sets poorly constrain spectral model parameters, and "best fit" models are often 
nonphysical. 

Investigators have adopted various approaches to analyze weak ACIS source spectra. Some use an 
assumed intrinsic source spectrum and absorbing column to derive a scaling factor between intrins ic lumi- 
nosity and observed flux that is applied to all weak sources |Hong et al.ll2005l : iMuno et alJl2006ah . Some 



2 CIAO calibrations include the effects of contamination on the ACIS-I optical blocking filter. 
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Fig. 3. — Light curves for some sources exhibiting strong variability during the Chandra observation, labeled 
by sequence number. Histograms show the full (0.5 — 8.0 keV) energy band with \/N error bars; dashed 
lines show the soft (0.5 — 2.0 keV) band; dotted lines show the hard (2.0 — 8.0 keV) band. Note that the 
ordinates vary between panels. 
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investigators apply prior knowledge ab out individual sources, such as Ay estimates from optical or NIR pho- 
tometry, to constrain spectral models ( Getman et al. 2006a ; Flaccomio et al. 20061) . Others limit their goals 
to obtaining non-unique spline-like fits to the event energy distribution, allowing nonphysical models, and 
emphasize caution when interpreting the spectral parameters (jFeigelson et al.ll2002l ) . All of these approaches 
are reasonable. 

Another issue concerns the mathematical method of seeking the best-fit spectral model, and then esti- 
mating confidence intervals on the spectral parameters. A least-squares approach based on the \ 2 statistic 
applied to binned data is traditional, although maximum likelihood estimation of unbinned data is becom- 
ing more common. With the nonlinear dependence between spectral parameters, it is often difficult to 
establish their confidence intervals for weak sources in an authoritative manner. Bayesian analysis pro- 
vides a likelihood-ba sed approach that incorporates prior knowledge expressed through prior distributions 



van Dvk et al.ll200lh . The uncertainties of the derived parameters (such as absorption-corrected luminosi- 



ties) are uncovered through posterior distributions, and one can marginalize over uninteresting parameters 
if one chooses. However, standard tools for a full Bayesian spectral analysis for X-ray data are not yet 
available. 



In this study, we adopt a strategy most similar to that of iFeigelson et al.l (|2002l ). Our primary goal 
is to obtain wide-band fluxes using regression models as spline functions. Fluxes are then converted to 
luminosities, both observed and intrinsic (corrected for estimated abso rption) , using a distance of 1.6 kpc. 
Calculations are performed using the XSP EC package (versio n 12.2.1o. lATnaulll99# l to fit one- and two- 
temperature apec thermal plasma models ([Smith et al.l 120011 ) and power law mo dels seen through an ab- 
sorbi ng column (Nh) of interstellar material with cosmic elemental abundances ( Morrison &: McCammonl 
1983). The one-temperature thermal plasma model is preferentially used and abundances of Z = O.SZp, a re 



assumed, as found for pre- main sequence stars in other, more nearby, star-forming regions (jlmanishi et al 



200 It IFeigelson et al.1 120021 ). A power law model is adopted if the thermal model poorly described the data 
or required nonphysical parameters, and if the source is not identified with a known stellar counterpart. 
When neither best-fit mo del is acceptable, we freeze the parameter kT = 2 keV in the thermal model, which 
is typical for young stars (jGetman et al. 2005b : Preibisch et al. 2005a ). and then fit for the Nh and normal- 
ization parameters. Analysis was not attempted on the 288 weakest sources whose photometric significance 
(column 12 in Tables CD and H]) is below 2.0. 

The b e st-fit model was foun d by the maximum likelihood method (XSPEC minimizes the C statistic, 
Cashlll979l IWachter et all 1 19791) . Although the accuracy of the parameter error estimation algorithm in 
XSPEC is not well understood for weak sources, we report 90% confidence intervals when they were flagged^ 
acceptable by XSPEC. One should regard the reported para meter errors with suitab le caution. A simulation 
of parameter uncertainties for weak sources was reported bv IFeigelson et al.1 (|2002l ). 



Spectral analysis results for the brighter 598 sources are presented in Tables EH and |H See the table 
notes for descriptions of the columns. Best-fit absorbing column densities range from negligible to log Nh ~ 
23.7 cm -2 , equivalent to a visual absorption of Ay ~ 250-300 mag. Temperatures range from the softest 
(kT ~ 0.4 keV) to the hardest (truncated at kT = 15 keV) accessible to the Chandra ACIS-I detectors. The 
range of observed total band (0.5 — 8 keV) luminosities, corrected for absorption, able to be derived from 
spectral modeling is 29.8 < logL t c < 33.3 ergs s _1 . Apparent luminosities range from \ogL t ~ 29.3 ergs s _1 



3 http : //heasarc . gsfc.nasa.gov/docs/ sof tware/lheasof t/xanadu/xspec 

4 See the telout error command in the XSPEC manual. 
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(estimated for the faintest "tentative" sources in Table [2]) to logLt = 32.8 ergs s 1 . Completeness limits are 
discussed in 



These tables report a variety of band-limited luminosities, both apparent and corrected for absorption 
(see Getman et al. 2005bl for details). Absorption-corrected soft-band luminosities (L s ,c m our nomen- 



clature) are highly uncertain and are not reported. For highly obscured sources (log Nh > 22.5 cm -2 ), 
soft-band luminosities are not reported because they may be misleading; the spectral fits are based solely on 
hard counts for these sources, so the model parameters may significantly underestimate or even miss entirely 
contributions from soft spectral components. Thus, although we could report L s and L StC derived from the 
model fit, we choose to omit these quantities because they are not representative of the true soft-band lumi- 
nosities of these obscured sources. The reader is cautioned that total-band luminosities may underestimate 
the true source luminosities due to such missing soft spectral components in our models. 

While a power law fit may suggest that the source could be a background AGN, this result is not 
definitive. Of the 52 sources in Table [4] all but 12 have IR counterparts (see Ej2.4|) ; brighter IR sources are 
less likely to be AGN. Conversely, some sources fit with very hard thermal plasmas may be equally well fit 
with power laws; we preferentially chose the thermal fit in such cases. 



2.4. Stellar Counterparts 

Due to the combination of high obscuration and confusion by bright nebular emission, there are no 
adequate visual band catalogs of stars in the M17 region. The best stellar survey of the region is in the 
NIR bands, using the SIRIUS instrument on the InfraRed Survey Facility (IRSF) 1.4 m telescope at the 
South African Astronomical Observatory. SIRIUS detected over 29,000 stars with s ensitivity limits a round 
J < 20 mag and K < 19 mag and sub-arcsecond resolution over a 14' x 14' region ( Jiang et al. 2002h . For 



outer portions of the ACIS field outside the SIRIUS coverage, we use the 2MASS All-Sky Survey which 
reaches K ~ 15.5 mag. A survey of the region with the IRAC detector on board the Spitzer Space Telescope 
in four IR (3.6 — 8 /im) bands has recently bee n released by the Galac tic Legacy Infrared Mid-Plane Survey 



Extraordinaire (GLIMPSE) Legacy Team (e.g. lChurchwell et al.ll2004f ). Here we considered the on-line ~ 30 



million source Highly Reliable Catalog from the GLIMPSE Enhanced Data Products@. 

We associate ACIS X-ray sources with IR sources using positional coincidence criteri^. We first removed 
reference frame offsets between the ACIS field (aligned to 2MASS) and the SIRIUS and GLIMPSE catalogs 
using high-quality ACIS sources (brighter than 10 counts and within 2' of the aim point). The estimated 
offsets were 0.6" to SIRIUS and 0.3" to GLIMPSE. Then, for every pair of catalog entries (e.g. an ACIS 
and a SIRIUS source), we test the hypothesis that the two physical sources are spatially coincident, i.e. that 
the two observed positions are random samples drawn from Gaussian distributions with identical means and 
with the reported standard deviations in position. If the hypothesis is true then the offsets between the two 
observed positions (in X and Y) will be drawn from a 2-D Gaussian distribution with zero mean and with 
X and Y variances equal to the sums of the reported variances. We reject this hypothesis if the observed 
positional offsets fall outside a 99% confidence region of this two-dimensional Gaussian. IR positional errors 
are obtained from the catalogs. ACIS positional errors consist of a random component (computed by AE), 



5 http : //www . astro .wise . edu/sirtf /glimpsedata. html 

6 Software implementing the matching algorithm described here is available in the TARA package at 
http : //www . astro . psu . edu/xray/docs/TARA/ . 
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corresponding to the uncertainty inherent in averaging a finite number of event positions, added in quadrature 
with an arbitrary 0.2" systematic component that avoids unreasonably small error estimates for very bright 
sources. 

The pairs of catalog entries for which the match hypothesis is satisfied are accepted as matches in order 
of their likelihood. When a match is accepted, the participating IR source is removed from all pending pairs 
not yet accepted, thus ensuring that no duplicates appear among the associations. In crowded locations, 
multiple IR sources satisfy the match hypothesis for a specific ACIS source. In such situations, if the most 
likely IR match of the group is still available, then it is accepted. Otherwise the ACIS source remains 
unmatched. 

Likely associations between ACIS sources and IR sources are reported in Table [5l Listed there are 
727 SIRIUS, 476 2MASS, and 224 GLIMPSE counterparts; 771 of the 886 ACIS sources (87%) have a 
counterpart identified. GLIMPSE counterparts were notably absent from the North Bar and South Bar 
which suffer from bright nebular emission. IR photometry is reproduced from the catalogs; 3 UK magnitudes 
are from SIRIUS photometry when available, except for four sources (ACIS #366, 701, 720, and 854) where 
2MASS values are more reliable due to s aturation effects in the SIRIUS im ages. Table [5] also gives likely 
associations to optical stars reported by Chini, Elsaesser. fc Neckel ( 1980l) (design ated CEN xxx) and to 
NIR stars from Hanson. Howarth. fc Contil ( 1997 ) derived from Bumgardner (1992) (designated B xxx) @. 
Association ambiguity flags and notes to other published characteristics of selected sources can be found in 
the table. 

The list of matches in Table [5] is not perfect, as both false negatives (true physical associations that are 
missed) and false positives (unphysical associations that are listed) are likely present. To the extent that the 
actual distributions of position errors follow the stated Gaussian distributions, we expect to miss 1% of true 
associations (false-negatives) as specified by the 99% confidence interval adopted. In Appendix[A]we estimate 
that 2.5% of the 2MASS associations, 7% of the SIRIUS associations, and 9% of the GLIMPSE associations 
are expected to be incorrect due to uncertainty in positions. A non-statistical problem in identifying stellar 
counterparts is that some cataloged IR sources may be condensations in the nebular emission (from atomic 
lines, PAH bands, or heated dust continuum) rather than stars. Some NIR sou rces may be unreso lved 
multiple systems, and (particularly when the primary has spectral type A or B, Istelzer et all l2005h the 
brightest component may not the principal X-ray emitter. In short, our associations should be considered 
only "likely" rather than confirmed stellar counterparts to the ACIS sources. 

Figure UJa shows the distributions of offsets between the ACIS and IR associations. Figure [4^ shows the 
expected deterioration in positional offsets as the Chandra PSF broadens off-axis. This can be compared 
to a similar plot from the Chandra Orion Ultradeep Project (COUP ) survey of the Orio n Nebula which 
shows median offsets increasing from 0.12" on-axis to 0.4" off-axis jGetman et al.ll2005bl Figure 9). The 
M17 median offsets are poorer, ranging from 0.2" to 1", due mainly to the much weaker ACIS sources in the 
short M17 observation compared to the long COUP observation. 



7 A reference frame offset of 1.9" between ACIS and lHanson. Howarth. fc Contil l ll997l) positions was removed before matching 
the catalogs. 
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Fig. 4. — Positional offsets between matching ACIS and IR sources, after removal of reference frame offsets. 
Panel (a) shows offset distributions for SIPJUS (dashed), 2MASS (solid), GLIMPSE (dotted), and Bum- 
gardncr (dash-dot) matches. The median ACIS/SIRIUS offset is 0.24". Panel (b) shows offsets for SIPJUS 
matches vs. off-axis angle in the ACIS field. 



Table 1. Primary Chandra Catalog: Basic Source Properties 



Source Position Extracted Counts Characteristics 



Seq 


C X O U J 


a J2000 


^J2000 














PSF 


Signif 


i°g P'b 






EffExp 


B median 


# 




(deg) 


(deg) 


K ) 


( ) 


Full 


Full 




Hard 


Frac 












(keV) 


f 1 1 














\P) 


\ 1 ) 


K°) 


{U) 


< i nl 


<ii) 












f 1 7'\ 


37 


182009.70-161503.7 


275 


040439 


-16 


251041 


0.8 


6.7 


9.7 


4.0 


2.3 


9.4 





90 


2.1 


<-5 




a 


32. 1 


5.3 


51 


182012.31-160447.7 


275 


051294 


-16 


079918 


0.8 


7.2 


16.4 


5.0 


3.6 


11.7 





90 


2.9 


<-5 




a 


31.6 


2.3 


97 


182017.84-161117.0 


275 


074374 


-16 


188070 


0.3 


3.1 


5.6 


3.0 


0.4 


5.7 





90 


1.6 


<-5 




a 


36.7 


3.2 


123 


182019.41-161329.1 


275 


080890 


-16 


224756 


0.1 


4.0 


30.8 


6.1 


0.2 


25.9 





60 


4.6 


<-5 




a 


36. 1 


3.3 


128 


182019.60-161038.8 


275 


081670 


-16 


177456 


0.3 


2.6 


8.6 


3.5 


0.4 


8.7 





89 


2. 1 


<-5 




a 


37.4 


4.8 


135 


182019.88-161052.3 


275 


082863 


-16 


181206 


0.3 


2.5 


8.5 


3.5 


0.5 


8.7 





89 


2.1 


<-5 




b 


37.4 


5.7 


163 


182020.73-160717.9 


275 


086416 


-16 


121661 


0.3 


3.9 


8.2 


3.5 


0.8 


7.5 





89 


2.0 


<-5 






36.4 


3.5 


177 


182021.21-161222.9 


275 


088397 


-16 


206382 


0.2 


2.9 


16.4 


4.6 


0.6 


12.6 





90 


3.1 


<-5 




a 


37.2 


2.8 


182 


182021.32-161140.2 


275 


088858 


-16 


194518 


0.2 


2.5 


11.5 


4.0 


0.5 


1.6 





89 


2.5 


<-5 




b 


37.4 


1.2 


186 


182021.43-160939. 1 


275 


089300 


-16 


160865 


0. 1 


2.3 


30.3 


6.1 


0.7 


16.6 





89 


4.6 


<-5 






38.0 


2.1 


187 


182021.43-161206.8 


275 


089325 


-16 


201912 


0.0 


2.7 


363.4 


19.6 


0.6 


288.6 





89 


18.1 


<-5 




a 


37.3 


3.1 


230 


182022.58-160251.3 


275 


094121 


-16 


047608 


0.2 


7.9 


365.4 


19.8 


5.6 


257.7 





90 


18.0 


<-5 




a 


31.2 


2.6 


233 


182022.69-160833.9 


275 


094578 


-16 


142761 


0.1 


2.7 


73.1 


9.1 


0.9 


24.6 





8 | ! 


7.6 


<-5 




a 


38.0 


1.7 


236 


182022.74-161123.5 


275 


094775 


-16 


189866 


0.1 


2.0 


60.3 


8.3 


0.7 


45.5 







6.8 


<-5 




a 


37.9 


2.6 


239 


182022.87-161148.5 


275 


095298 


-16 


196825 


0.1 


2.2 


19.3 


5.0 


0.7 


19.5 





89 


3.5 


<-5 




a 


37.1 


4.9 


244 


182022.93-161152.7 


275 


095573 


-16 


197986 


0.1 


2.3 


134.4 


12.1 


0.6 


77.6 





89 


10.6 


<-5 






35. 1 


2.4 


246 


182022.97-161131.8 


275 


095712 


-16 


192185 


0.2 


2.1 


12.4 


4.1 


0.6 


12.6 





^ ' 


2.6 


<-5 


g 




34.1 


4.8 


255 


182023.16-161305.7 


275 


096504 


-16 


218256 


0.1 


3.1 


54.4 


7.9 


0.6 


54.5 





90 


6.4 


<-5 




a 


35.5 


4.6 


275 


182023.80-161325.7 


275 


099198 


-16 


223807 


0.4 


3.3 


8.5 


3.5 


0.5 


8.6 





90 


2.1 


<-5 




a 


34.8 


4.9 


281 


182024.00-160818.8 


275 


100040 


-16 


138573 


0.1 


2.7 


29.2 


6.0 


0.8 


1.7 





89 


4.5 


<-5 




a 


36.2 


0.9 


296 


182024.39-160843.3 


275 


101640 


-16 


145367 


0.2 


2.3 


10.2 


3.8 


0.8 


1.7 





90 


2.3 


<-5 






38.4 


1.3 


309 


182024.60-161139.2 


275 


102519 


-16 


194249 


0.1 


1.8 


19.0 


5.0 


1.0 


18.3 





90 


3.4 


<-5 




a 


38.2 


3.3 


324 


182024.87- 161127.5 


275 


103634 


-16 


190973 


0. 1 


1.6 


49.8 


7.7 


1.2 


31.2 





90 


6.1 


<-5 




a 


38.2 


2.7 


329 


182024.94-161131.1 


275 


103949 


-16 


191974 


0.2 


1.7 


15.7 


4.6 


1.3 


16.2 





90 


3.0 


<-5 




c 


38.2 


4.4 


334 


182025.07—161133.9 


275 


104499 


-16 


192767 


0.1 


1.7 


21.8 


5.3 


1.2 


12.2 





90 


3.7 


<-5 




a 


38.2 


2.5 


350 


182025.50-161053.8 


275 


106253 


-16 


181624 


0.0 


1.2 


308.6 


18.1 


1.4 


188.5 





90 


16.6 


<-5 






38.6 


2.5 


370 


182025.93-161114.3 


275 


108062 


-16 


187315 


0.2 


1.3 


8.8 


3.7 


1.2 


1.6 





90 


2.1 


<-5 




a 


38.6 


1.5 


371 


182025.93-160949.4 


275 


108070 


-16 


163734 


0.2 


1.3 


8.1 


3.5 


0.9 


3.6 





90 


2.0 


<-5 




a 


38.9 


1.8 


372 


182025.95-160938.2 


275 


108153 


-16 


160626 


0.3 


1.4 


5.1 


3.0 


0.9 


3.6 





90 


1.4 


-3.6 




b 


38.8 


2.5 


373 


182026.00-161253.0 


275 


108352 


-16 


214727 


0.3 


2.6 


5.6 


3.0 


0.4 


5.8 





85 


1.6 


<-5 






37.8 


4.0 


374 


182026.01-161022.4 


275 


108377 


-16 


172897 


0.0 


1.1 


102.3 


10.7 


0.7 


46.7 





90 


9.1 


<-5 






27.4 


1.8 


375 


182026.02-161240.0 


275 


108417 


-16 


211132 


0.2 


2.4 


8.4 


3.5 


0.6 


4.6 





89 


2.0 


<-5 




a 


38.0 


2.3 


376 


182026.04-161104.6 


275 


108512 


-16 


184620 


0.1 


1.2 


78.3 


9.5 


1.7 


45.3 





90 


7.8 


<-5 




a 


38.6 


2.3 


377 


182026.09-161039.9 


275 


108747 


-16 


177772 


0.2 


1.0 


13.9 


4.4 


1.1 







90 


2.8 


<-5 




a 


38.7 


1.8 


378 


182026.09-160301.3 


275 


108749 


-16 


050381 


0.7 


7.5 


11.5 


4.5 


4.5 


0.1 





90 


2.3 


-4.7 




a 


31.8 


1.1 


379 


182026. 11- 161053.5 


275 


108806 


-16 


181547 


0.1 


1.1 


21.0 


5.3 


2.0 


10.1 





90 


3.6 


<-5 




a 


38.7 


1.9 


396 


182026.60-161055.7 


275 


110852 


-16 


182141 


0.0 


1.0 


2473.5 


50.3 


2.5 


1622.0 


° 




48.7 


<-5 




c 


38.7 


2.7 


398 


182026.62-161136.5 


275 


110920 


-16 


193480 


0.1 


1.4 


18.2 


4.9 


0.8 


7.5 




90 


3.3 


<-5 




a 


38.4 


1.9 


399 


182026.62-160822.9 


275 


110937 


-16 


139700 


0.3 


2.3 


6.2 


3.2 


0.8 


1.7 





90 


1.6 


-4.8 






36.0 


1.7 


433 


182027.41-161331.0 


275 


114226 


-16 


225280 


0.0 


3.1 


263.3 


16.8 


0.7 


18.6 





90 


15.2 


<-5 




a 


37.5 


1.3 


466 


182028.15-161049.3 


275 


117308 


-16 


180378 


0.0 


0.6 


538.1 


23.7 


1.9 


386.2 





90 


22.2 


<-5 




a 


38.9 


2.7 


488 


182028.65-161211.6 


275 


119414 


-16 


203233 


0.1 


1.7 


70.6 


8.9 


0.4 


35.8 





89 


7.5 


<-5 




a 


38.5 


2.0 


495 


182028.86-161042.0 


275 


120280 


-16 


178337 


0.2 


0.4 


8.3 


3.9 


2.7 


0.0 





90 


1.9 


-4.0 






39.0 


1.6 


510 


182029.21-161425.6 


275 


121709 


-16 


240461 


0.4 


3.9 


5.2 


3.0 


0.8 


0.5 





90 


1.4 


-3.8 






36.7 


1.2 


511 


182029.28-161041.4 


275 


122029 


-16 


178180 


0.1 


0.3 


53.0 


8.2 


6.0 


29.7 





90 


6.1 


<-5 






39.1 


2.2 


512 


182029.31-161047.9 


275 


122151 


-16 


179989 


0.1 


0.4 


77.4 


9.5 


3.6 


50.4 





83 


7.7 


<-5 




a 


39.1 


2.5 


513 


182029.32-161045.2 


275 


122167 


-16 


179236 


0.1 


0.4 


16.8 


5.0 


3.2 


11.5 





79 


3.0 


<-5 




a 


39.1 


2.4 


514 


182029.36-160919.2 


275 


122353 


-16 


155335 


0.2 


1.2 


6.1 


3.2 


0.9 


2.6 





91 


1.6 


-4.4 




a 


39.1 


1.8 


515 


182029.39-160943.0 


275 


122474 


-16 


161946 


0.0 


0.8 


155.2 


13.0 


0.8 


102.6 





90 


11.5 


<-5 




a 


39.3 


2.5 


516 


182029.39-161046.5 


275 


122489 


-16 


179593 


0.0 


0.4 


69.5 


9.1 


3.5 


53.4 





79 


7.2 


<-5 




a 


39.1 


3.0 


517 


182029.39-161212.9 


275 


122489 


-16 


203593 


0.2 


1.7 


8.5 


3.5 


0.5 


1.8 





89 


2.1 


<-5 




a 


38.5 


1.4 


518 


182029.43-161050.2 


275 


122661 


-16 


180616 


0.1 


0.4 


47.0 


7.8 


6.0 


29.9 





90 


5.6 


<-5 




a 


39.1 


2.7 


536 


182029.81-161045.6 


275 


124213 


-16 


179352 


0.0 


0.3 


1912.2 


44.3 


6.8 


1052.5 





80 


42.7 


<-5 




a 


39.1 


2.2 


543 


182029.89-161044.5 


275 


124574 


-16 


179031 


0.0 


0.3 


3865.9 


62.7 


8.1 


2447.0 





80 


61.1 


<-5 






39.1 


2.6 


567 


182030.23-161034.9 


275 


125982 


-16 


176382 


0.1 


0.1 


43.6 


7.5 


4.4 


9.4 





90 


5.5 


<-5 






39.2 


1.7 


574 


182030.44-161053.1 


275 


126836 


-16 


181427 


0.1 


0.4 


89.5 


10.3 


5.5 


22.3 





90 


8.3 


<-5 






39.2 


1.7 


600 


182030.95-161039.4 


275 


128970 


-16 


177620 


0.2 


0.2 


8.0 


3.7 


2.0 


5.0 





90 


1.9 


-4.3 




a 


39.2 


2.3 


618 


182031.35-160228.4 


275 


130660 


-16 


041249 


0.4 


8.0 


64.4 


9.3 


11.6 


38.7 





75 


6.5 


<-5 




a 


29.3 


2.3 



Table 1 — Continued 





Source 




Position 








Extra 


otcd Coi 


ints 








Characteristics 






Scq 


CXOU J 


°J2000 


d .J2000 


Err 


e 


Net 


ANct 


Bkgd 


Net 


PSF 


Sigrrif 


logP B 


Arrom Var 


Eff'Exp 




# 




(deg) 


(deg) 


(") 


(') 


Pull 


Full 


Full 


Hard 


Frac 








(ks) 


(keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) (15) 


(16) 


(17) 


647 


182031.84-161138.1 


275.132701 


-16.193943 


0.3 


1.2 


3.5 


2.5 


0.5 


0.8 


0.90 


1.1 


-2.9 


g 


27.5 


1.6 


650 


182031.89-161616.5 


275.132883 


-16.271262 


0.6 


5.8 


10.0 


3.8 


1.0 


3.5 


0.91 


2.3 


<-5 


g 


18.0 


1.7 


655 


182031.97-161030.5 


275.133240 


-16.175147 


0.0 


0.4 


249.3 


16.3 


0.7 


208.6 


0.90 


14.8 


<-5 


g 


25.8 


3.4 


706 


182034.56-161523.6 


275.144010 


-16.256565 


0.5 


5.0 


9.0 


3.7 


1.0 


9.4 


0.89 


2.1 


<-5 




33.6 


5.0 


726 


182035.63-161055.5 


275.148459 


-16.182086 


0.2 


1.3 


11.2 


4.0 


0.8 


1.8 


0.90 


2.5 


<-5 




39.0 


1.4 


854 


182053.86-160306.5 


275.224439 


-16.051830 


0.4 


9.3 


105.5 


11.3 


9.5 


6.0 


0.91 


8.9 


<-5 




29.3 


1.0 



Note. — Tabic is published in its entirety in the electronic edition of the Astrophysical Journal. Interesting sources mentioned in the text are shown here for guidance regarding its form and 

Note. — Column 1: X-ray catalog sequence number, sorted by RA. Column 2: IAU designation. Columns 3,4: Right ascension and declination for epoch J2000.0. Column 5: Estimated 
standard deviation of the random component of the position error, ^/tr^ + <x^ . The single-axis position errors, cr x and cry . arc estimated from the single-axis standard deviations of the PSF inside 
the extraction region and the number of counts extracted. Column 6: Off-axis angle. Columns T,8: Net counts extracted in the total energy band (0.5—8 keV); average of the upper and lower 
Ict errors on column 7. Column 9: Background counts extracted (total band). Column 10: Net counts extracted in the hard energy band (2-8 keV). Column 11: Fraction of the PSF (at 1.497 
keV) enclosed within the extraction region. Note that a reduced PSF fraction (significantly below 90%) may indicate that the source is in a crowded region. Column 12: Photometric significance 
computed as uppcr c ™^r C °n net counts ' Column 13: Lo S probability that extracted counts (total band) arc solely from background. Some sources have P B values above the 1% threshold that 
defines the catalog because local background estimates can rise during the final extraction iteration after sources are removed from the catalog. Column 14: Source anomalies: g = fractional time 
that source was on a detector (FRACEXPO from mkarf) is < 0.9 ; c = source on field edge; p = source piled up; s = source on readout streak. Column 15: Variability characterization based on 
K-S statistic (total band): a = no evidence for variability (0.05 < PkS^' b = P oss ibly variable (0.005 < P K g < 0.05); c = definitely variable (Pj^s < 0.005). No value is reported for sources with 
fewer than 4 counts or for sources in chip gaps or on field edges. Column 16: Effective exposure time: approximate time the source would have to be observed on-axis (no telescope vignetting) on 
a nominal region of the detector (no dithering over insensitive regions of the detector) to obtain the reported number of counts. Column 17: Background-corrected median photon energy (total 



Table 2. Tentative Chandra Catalog: Basic Source Properties 





Source 




Position 








Extracted Counts 








Characteristics 






Seq 


CXOU J 


Q J2000 


5 J2000 


Err 





Net 


ANct 


Bkgd 


Net 


PSF 




log P B 


Atom Var 


EffExp 


E median 


# 




(deg) 


(deg) 


(") 


o 


Full 


Full 


Full 


Hard 


Frac 








(ks) 


(kcV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(') 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) (15) 


(16) 


(17) 


322 


182024.83-161135.3 


275.103482 


-16.193161 


0.5 


1.7 


2.9 


2.5 


1.1 


3.3 


0.90 


0.9 


-1.6 




38.2 


4.7 


519 


182029.48-161644.0 


275.122866 


-16.278895 


0.7 


6.2 


5.6 


3.4 


2.4 


0.6 


0.90 


1.4 


-2.5 




33.7 


1.5 


Note. 


— Table [2] has the san 


le columns as ^] and 


is published 




entirety i 


n the electronic edition of the 


A Htrophystcal Jo u.nial. 


Intere 




s mentioned in the 




shown here for 



guidance regarding its form and content. 



Note. — Sec notes for Tablc^ 



Table 3. X-ray Spectroscopy for Photometrically Selected Sources: Thermal Plasma Fits 



Sourcc a 



Spectral Fit b 



X-ray Luminosities 



Scq 


CXOU J 




Not 


Signif 


1 








kT 




k 


tgEM 




logL s 


log 




logife.c 


log 




log it, 


# 








Counts 




(cm 


-2 


) 






(koV) 




(cm-" 


) 












(ergs s -1 ) 








(1) 




(2) 




(3) 


(4) 




(5) 








(6) 






(7) 






(8) 


0) 


(10) 


(11) 


(12) 


37 


182009 


70-161503 


7 


9.7 


2.1 




23 


7 








1.4 






55.7 










30 


55 


31 


99 


30 


55 




51 


182012 


31-160447 


7 


16.4 


2.9 


-0.3 


22 


3 


+0 


2 




2.0 






53.9 


+0 


2 


29 


54 


30 


27 


30 


41 


30 


35 


30.85 


123 


182019 


41-161329 


1 


30.8 


4.6 


-0.3 


22 


5 


+ 


2 


-1.7 3.4 


-0 


4 


54.2 


+0 


5 


29 


68 


30 


86 


31 


01 


30 


89 


31.30 


128 


182019 


60-161038 


8 


8.6 


2.1 


-0.5 


23 


2 


+ 


3 




13.6 






53.8 










30 


45 


30 


82 


30 


45 




135 


182019 


88-161052 


3 


8.5 


2.1 




23 


8 








2.0 






55.3 










30 


54 


31 


85 


30 


54 




177 


182021 


21-161222 


9 


16.4 


3.1 


-0.6 


22 


3 


+ 


4 


— 2 


3 3.6 






53.6 


+0 


7 


29 


40 


30 


32 


30 


42 


30 


37 


30.71 


182 


182021 


32-161140 


2 


11.5 


2.5 


-1.1 


21 


9 


+ 


3 


-0 


4 0.7 +0.9 






53.5 


+0 


9 


29 


66 


29 


06 


29 


14 


29 


76 


30.47 


186 


182021 


43-160939 


1 


30.3 


4.6 


-0.6 


22 


1 


+ 


4 


-1.7 3.0 


-0 


4 


53.8 


+0 


6 


29 


81 


30 


42 


30 


48 


30 


52 


30.81 


187 


182021 


43-161206 


8 


363.4 


18.1 


-0.08 


22 


4 


+ 


07 


-l 


2 4.4 +2.6 


-0 


1 


55.0 


+0 


1 


30 


56 


31 


76 


31 


89 


31 


79 


32.14 


230 


182022 


58-160251 


3 


365.4 


18.0 


-0.10 


22 


3 


+0 


09 


-l 


5 4.6 +3.3 


-0 


1 


55.0 


+0 


1 


30 


77 


31 


74 


31 


83 


31 


79 


32.08 


233 


182022 


69-160833 


9 


73.1 


7.6 


-0.4 


21 


7 


+0 


2 


-l. 


5 3.5 +10.7 


-0 


2 


53.9 


+0 


2 


30 


29 


30 


69 


30 


72 


30 


84 


31.01 


236 


182022 


74-161123 


5 


60.3 


6.8 


-0.2 


22 


6 


+0 


2 


-0 


5 1.6 +1.8 


-0 


5 


54.7 


+0 


4 






30 


86 


31 


12 


30 


90 




239 


182022 


87-161148 


5 


19.3 


3.5 


-0.2 


23 


5 


+ 


2 




2.0 


-0 


3 


55.2 


+0 


4 






30 


78 


31 


69 


30 


78 




244 


182022 


93-161152 


7 


134.4 


10.6 


-0.2 


22 


2 


+ 


1 


-1 


1 3.3 +7.5 


-0 


3 


54.6 


+0 


2 


30 


42 


31 


24 


31 


33 


31 


30 


31.63 


246 


182022 


97-161131 


8 


12.4 


2.6 


-0.3 


23 


5 


+ 


3 




2.6 






54.9 


+1 


6 






30 


72 


31 


56 


30 


72 




255 


182023 


16-161305 


7 


54.4 


6.4 


-0.3 


23 


5 


+ 


2 


-0 


7 1.6 +8.5 






55.8 










31 


20 


32 


18 


31 


20 




275 


182023 


80-161325 


7 


8.5 


2.1 


-0.4 


23 


1 


+ 


3 




13.7 






53.8 










30 


44 


30 


79 


30 


44 




281 


182024 


00-160818 


8 


29.2 


4.5 




20.0 






-0 


3 0.6 +o.i 


-0 


2 


53.3 


+0 


2 


30 


21 


28 


80 


28 


80 


30 


22 


30.24 


296 


182024 


39-160843 


3 


10.2 


2.3 




21 


6 


+0 


4 




2.0 


-0 


4 


53.1 


+0 


3 


29 


52 


29 


58 


29 


60 


29 


85 


30.04 


309 


182024 


60-161139 


2 


19.0 


3.4 


-0.3 


22 


8 


+ 


4 


-1 


1 1.7 +3.4 






54.5 


+1 


7 






30 


50 


30 


88 


30 


51 




324 


182024 


87-161127 


5 


49.8 


6.1 


-0.2 


22 


4 


+0 


2 


-1 


o 2.5 +4.o 


-0 


3 


54.2 


+0 


4 


29 


88 


30 


74 


30 


87 


30 


80 


31.24 


329 


182024 


94-161131 


1 


15.7 


3.0 




23 


1 








15.0 






54.0 










30 


66 


31 


00 


30 


66 




334 


182025 


07-161133 


9 


21.8 


3.7 


-0.3 


22 


4 


+ 


3 


-0 


8 1.7 +2.9 






54.1 






29 


61 


30 


30 


30 


47 


30 


38 


30.99 


350 


182025 


50-161053 


8 


308.6 


16.6 


-0.12 


22 


1 


+0 


10 


-1 


6 5.0 +5.6 


-0 


1 


54.7 


+0 


1 


30 


72 


31 


55 


31 


62 


31 


61 


31.85 


370 


182025 


93-161114 


3 


8.8 


2.1 




21 


8 








2.4 






53.0 






29 


28 


29 


60 


29 


64 


29 


77 


30.03 


374 


182026 


01-161022 


4 


102.3 


9.1 


-0.3 


22 





+0 


2 


-1 


3 3.3 +7.3 


-0 


2 


54.4 


+0 


2 


30 


49 


31 


09 


31 


14 


31 


19 


31.45 


375 


182026 


02-161240 





8.4 


2.0 


-0.6 


22 


6 


+ 


3 


-0 


7 1.3 






54.1 


+1 


4 






29 


98 


30 


29 


30 


03 




376 


182026 


04-161104 


6 


78.3 


7.8 


-0.2 


22 


3 


+0 


3 


-1 


1 1.9 +1.6 


-0 


3 


54.4 


+0 


7 


30 


20 


30 


81 


30 


93 


30 


91 


31.38 


377 


182026 


09-161039 


9 


13.9 


2.8 


-0.6 


22 





+ 


5 


-1 


6 2.2 






53.4 


+1 


2 


29 


51 


29 


91 


29 


97 


30 


06 


30.39 


378 


182026 


09-160301 


3 


11.5 


2.3 




20 


4 


+1 


6 




1.4 


-0 


3 


53.1 


+0 




29 


84 


29 


36 


29 


36 


29 


96 


29.98 


379 


182026 


11-161053 


5 


21.0 


3.6 


-0.4 


22 





+ 


3 


-0 


7 1.8 +3.1 


-0 


4 


53.7 


+ 


4 


29 


73 


30 


07 


30 


14 


30 


23 


30.62 


396 


182026 


60-161055 


7 


2473.5 


48.7 


-0.04 


22 


1 


+0 


03 


-4. 


i 14.5 


-0.02 


55.6 


+0 


02 


31 


56 


32 


54 


32 


59 


32 


58 


32.75 


398 


182026 


62-161136 


5 


18.2 


3.3 


-0.2 


22 


5 


+0 


2 


-0 


3 0.7 +0.5 






54.6 


+1 









29 


99 


30 


30 


30 


17 




433 


182027 


41-161331 





263.3 


15.2 


-0.07 


22 





+0 


08 


-0.11 0.6 +0.08 


-0 


2 


55.1 


+0 


3 


31 


02 


30 


32 


30 


43 


31 


10 


32.06 


466 


182028 


15-161049 


3 


538.1 


22.2 


-0.10 


22 


3 


+ 


07 


-1 


4 4.8 +5.2 


-0 


1 


55.1 


+0 


1 


30 


85 


31 


86 


31 


95 


31 


91 


32.19 


488 


182028 


65-161211 


6 


70.6 


7.5 


-0.3 


22 


2 


+ 


2 


-0 


8 2.2 +2.7 


-0 


3 


54.3 


+0 


3 


30 


15 


30 


78 


30 


89 


30 


87 


31.30 


511 


182029 


28-161041 


4 


53.0 


6.1 




22 


6 








0.9 +0.5 


-0 


6 


55.0 










30 


55 


30 


89 


30 


65 




512 


182029 


31-161047 


9 


77.4 


7.7 


-0.2 


22 


5 


+ 


2 


-0 


5 1.5 +1.0 


-0.4 


54.8 


+ 


5 






30 


95 


31 


19 


31 


00 




513 


182029 


32-161045 


2 


16.8 


3.0 


-0.2 


22 


5 


+0 


2 




2.0 


-0 


3 


54.0 


+0 


2 


29 


43 


30 


34 


30 


52 


30 


39 


30.96 



Notes 



(13) 



Table 3 — Continued 



Source a Spectral Fit b X-ray Luminosities Notes 



Scq 


CXOU J 




Net 


Signif 


log N H 


kT 


log EM 


log 




log 






log 




log it, c 




# 






Counts 




(cm" 2 


) 


(keV) 


(cm- 3 


) 










(ergs s- 1 ) 












(1) 


(2) 




(3) 


(4) 


(5) 




(6) 


(7) 




(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


515 


182029.39-160943 





155.2 


11.5 


-0.2 22.3 


+ 0.1 


-0.8 2.5 +2.9 


-0.3 54.7 




30 


39 


31 


23 


31 


36 


31 


29 


31 


72 




516 


182029.39-161046 


5 


69.5 


7.2 


-0.2 22.5 


+ 0.2 


-1.5 3.0 + 7.7 


-0.3 54.5 


+ 0.5 


29 


95 


31 


04 


31 


18 


31 


07 


31 


51 




517 


182029.39-161212 


9 


8.5 


2.1 


-0.7 22.2 


+0.3 


0.4 


54.3 




29 


48 


28 


94 


29 


14 


29 


59 


31 


13 




518 


182029.43-161050 


2 


47.0 


5.6 


22.5 




1.5 


54.5 


+0.1 






30 


66 


30 


90 


30 


72 








536 


182029.81-161045 


6 


1912.2 


42.7 


-0.03 22.4 


+ 0.03 


15.0 


-0.04 55.3 


+ 0.04 


31 


63 


32 


37 


32 


48 


32 


44 


33 


16 


2T (0.9 kcV) 


543 


182029.89-161044 


5 


3865.9 


61.1 


-0.05 22.3 


+ 0.04 


-1.9 10.4 +3.1 


-0.02 55.7 


+ 0.02 


31 


85 


32 


78 


32 


85 


32 


83 


33 


28 


2T (0.7 kcV) 


567 


182030.23-161034 


9 


43.6 


5.5 


-o.i 22.3 


+0.2 


-0.3 0.6 +0.3 


-0.4 54.8 


+ 1.2 


30 


18 


30 


05 


30 


26 


30 


42 


31 


72 




574 


182030.44-161053 


1 


89.5 


8.3 


22.4 




0.6 +0.2 


-0.4 55.3 




30 


47 


30 


42 


30 


69 


30 


75 


32 


25 




618 


182031.35-160228 


4 


64.4 


6.5 


22.5 




1.3 


-0.59 54.9 


+ 0.08 


30 


25 


30 


86 


31 


09 


30 


96 


31 


81 




650 


182031.89-161616 


5 


10.0 


2.3 


20.6 




14.9 


53.0 




29 


69 


30 


07 


30 


08 


30 


23 


30 


24 




655 


182031.97-161030 


5 


249.3 


14.8 


-0.14 22.5 


+ 0.09 


-4.0 8.9 


-0.08 54.9 


+ 0.13 


30 


46 


31 


83 


31 


95 


31 


85 


32 


13 


1 


706 


182034.56-161523 


6 


9.0 


2.1 


-0.6 23.3 


+0.3 


15.0 


-o.e 54.0 


+ 0.4 






30 


56 


31 


02 


30 


56 








726 


182035.63-161055 


5 


11.2 


2.5 


-0.8 21.8 


+ 0.3 


2.0 


-0.3 53.2 


+ 0.3 


29 


51 


29 


71 


29 


75 


29 


92 


30 


19 


~J 


854 


182053.86-160306 


5 


105.5 


8.9 


21.0 


+ 0.2 


0.6 +o.i 


54.0 




30 


83 


29 


45 


29 


47 


30 


85 


30 


99 


1 



Note. — Table [3] is published in its entirety in the electronic edition of the Astrophysical Journal. Interesting sources mentioned in the text are shown here for guidance regarding 
its form and content. 

a For convenience, columns 1—4 reproduce the source identification, net counts, and photometric significance data from Table IT1 

b All fits used the "wabs(apec)" model in XSPEC and assumed 0.3.Zq abundances Jlmanishi et alll2Q0ll; )Feigelson et afl|2002|) . Columns 5 and 6 present the best-fit values for the 
column density and plasma temperature parameters. Column 7 presents the emission measure for the model spectrum, assuming a distance of 1.6 kpc. Quantities in italics were 
frozen in the fit. Uncertainties represent 90% confidence intervals. More significant digits are used for uncertainties < 0.1 in order to avoid large rounding errors; for consistency, the 
same number of significant digits is used for both lower and upper uncertainties. Uncertainties arc missing when XSPEC was unable to compute them or when their values were so 
large that the parameter is effectively unconstrained. Fits lacking uncertainties, fits with large uncertainties, and fits with frozen parameters should be viewed merely as splines to the 
data to obtain rough estimates of luminosities; the listed parameter values are unreliable. 

c X-ray luminosities arc presented in columns 8—12: s — soft band (0.5-2 keV); h — hard band (2—8 keV); t — total band (0.5—8 keV). Absorption-corrected luminosities arc 
subscripted with a c; they are omitted when log Nh > 22.5 cm" since the soft band emission is essentially unmeasurable. 

d 2T means a two-temperature model was used; the second temperature is shown in parentheses. See spectra in Figure [121 



Table 4. X-ray Spectroscopy for Photometrically Selected Sources: Power Law Fits 





Source 3 








Spectral Fit b 














X-ray 


Fluxes 








Soq 


CXOU J 




Not 


Signif 


log N H 


r 


log 


N r 


log 


L s 


log 






log 


Lt 


log L t ,c 


# 






Counts 




(cm" 2 ) 
















(photons 


cm 2 s 


— 1- 


1 




(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


2 


181957.67-161328. 


.8 


12.5 


2.3 


23.1 


1. 


1 


-6 


,3 






29 


.82 


30 


.10 


29 


.82 




5 


181959.04-160350 


2 


55.4 


5.7 


21.9 


1. 


6 


-5 


.2 


30. 


10 


.30 


.88 


30 


91 


.30 


.95 


31.10 


13 


182004.28-161128. 


9 


19.8 


3.4 


21.6 


-1.3 2. 


.1 


-5 


.6 


29 


.78 


30 


.12 


30 


.14 


.30 


.28 


30.46 


40 


182010.53-161401. 


.5 


24.1 


3.9 


-l.o 21.9 +0.4 


-1.0 1. 


.5 


-5 


.6 


29. 


70 


30 


.49 


30 


.52 


30 


.50 


30.71 


59 


182013.22-160856 


.5 


10.3 


2.3 


22.2 


0. 


1 


-6 


.4 


28. 


95 


30 


.04 


30 


.67 


30 


.65 


30.70 


69 


182015.48-160706. 





10.1 


2.3 


22.8 


5. 


.4 


-3 


.6 






29 


.97 


30 


.54 


30 


.02 




96 


182017.83-161014. 


.5 


12.8 


2.7 


-1.4 22.8 +0.5 


1. 


2 


-5 


.3 






30 


.90 


31 


.08 


30 


.90 




102 


182018.19-160937. 


.0 


56.5 


6.6 


-0.3 22.5 +0.2 


-1.0 1. 


.5 


-4 


.9 






31 


.07 


.31 


.18 


31 


.08 




112 


182018.95-160736. 


9 


17.3 


3.2 


-0.7 22.3 +0.4 


1. 


9 


-5 


.4 


29 


33 


30 


.41 


30 


.50 


.30 


.45 


30.76 


121 


182019.33-161051. 


9 


9.5 


2.2 


23.0 


3 


.1 


-4 


.4 






30 


32 


30 


.78 


30 


.33 




153 


182020.49-160806 


4 


15.0 


2.9 


22.0 +0.5 


1. 


6 


-5 


.7 


29 


45 


30 


.30 


30 


.34 


30 


.35 


30.54 


175 


182021.08-160738. 


2 


20.2 


3.6 


22.0 +0.5 


1. 


3 


-5 


.8 


29 


47 


30 


.50 


30 


.54 


30 


.54 


30.67 


179 


182021.25-160944. 


3 


8.4 


2.0 


21.5 +0.8 





6 


-6 


.6 


29 


11 


30 


.18 


30 


.19 


.30 


.21 


30.24 


189 


182021.59-160703. 


.3 


11.3 


2.5 


22.5 


1. 


6 


-5 


.6 






30 


.30 


30 


.43 


30 


.32 




196 


182021.74-161317. 


4 


10.4 


2.3 


22.0 +0.5 


1. 


.5 


-6 


.0 


29. 


24 


30 


.16 


30 


.20 


30 


.21 


30.37 



Note. — Tabled is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding 
its form and content. 

a For convenience, columns 1—4 reproduce the source identification, net counts, and photometric significance data from Table |T1 

b All fits used the "wabs(powcrlaw)" model in XSPEC. Columns 5 and 6 present the best-fit values for the column density and power law photon 
index parameters. Column 7 presents the power law normalization for the model spectrum. Quantities in italics were frozen in the fit. Uncertainties 
represent 90% confidence intervals. Uncertainties are missing when XSPEC was unable to compute them or when their values were so large that the 
parameter is effectively unconstrained. Fits lacking uncertainties, fits with large uncertainties, and fits with frozen parameters should be viewed merely 
as splines to the data to obtain rough estimates of luminosities; the listed parameter values arc unreliable. 



c X-ray luminosities arc presented in columns 8—12: s — soft band (0.5-2 kcV); h — hard band (2-8 keV); t — total band (0.5—8 keV). Absorption- 
corrected luminosities arc subscripted with a c; they are omitted when log Nh > 22.5 cm -2 since the soft band emission is essentially unmcasurablc. 



Table 5. Stellar counterparts 





X-ray 








Optical/I 


nfrar 


ed star 










Infra 


cd photor 


letry 


Scq 


cxou 




CEN 


B 


2MASS 




SIRIUS 


GLIMPSE 


Notes 


J 

(mag) 


H 

(mag) 


K 

(mag) 


[3.6] 

(mag) 


[4-5] 

(mag) 


(1) 


(2) 




(3) 


(4) 


(5) 




(6) 


(7) 


(S) 


(9) 


(10) 


(11) 


(12) 


(13) 


37 


182009.70 — 161503 


7 


























51 


182012.31 — 160447 


7 






18201227—1604466 


Obc 


25643 


G015.1 109 — 00.5795 


051 


15.89 


12.27 


10.35 


8.70 


8.39 


97 


182017.84—161117 













13045 




097 




17.81: 


16.47: 






123 


182019.41 — 161329 


1 










15128 




123 


12.75 


9.68 


7.45 






128 


182019.60-161038 


8 




353 


18201960—1610382 
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2.4-1- Sources Without Stellar Counterparts 

IR associations could not be found for 115 of the 886 X-ray sources. Their spatial locations are plotted 
as small diamonds in Figure [TJa. 

Visual examination of the SIRIUS and 2MASS images suggests that at least 14 of these probably have 
NIR associations that are not cataloged, either because the counterpart is only marginally resolved from a 
brighter star or because it lies near the limit of the NIR survey. Visual examination of the GLIMPSE image 
shows many point-like objects that are not listed in the Highly Reliable Catalog; thus some ACIS sources 
unidentified in Table [5] may have 3.6 — 8 fim counterparts. 

Extragalactic sources, mostly AGN at moderate redshifts, dominate X-ray source counts a t high Galactic 



latitud es an d some shoul d be de tected in Galactic plane fields despite the heavy obscuration. iGetman et al 



(|2006af ) and IWang et al.l (|2006f ) have quantitatively calculated the expected background contamination in 
ACIS observations of the Cepheus B and NGC 6357 star formation regions. In both cases, they predict that 
20 ± 10 sources should appear in a ~40 ks ACIS exposure with dense Galactic cloud material, and only a few 
of those should be detectable in the IR. Some of these extragalactic sources are probably among the sources 
with spectra best fit by power law rather than thermal models (Table 2]); as noted above, twelve of those 
sources lack IR counterparts. The spatial distribution of unassociated sources in Figure Q]a is consistent with 
a random sprinkling of ~20 extragalactic sources. 

The remaining ~80 ACIS sources without counterparts are likely to be new members of the M17 cloud 
complex. The high median energy of this unassociated population (E mec n an — 3.0 keV) indicates that these 
sources are deeply embedded (Ay ^15 mag). Some are probably low mass stars with strong magnetic 
activity; several d ozen such X-ray disco vered stars were found in the cloud behind the Orion Nebula in 



the COUP study (jGetman et alJ l2005af ). This interpretation is supported by their spatial distribution: 
Figure Q]a shows that they are concentrated in the IR-bright North Bar and South Bar where IR surveys 
are most confused by nebular emission. Some of these may be very young with the local absorption in an 
envelope or disk, as suggested by Getman et al. for the COUP sources. 



24.2. NIR Properties of ACIS Sources 



Following Hanson. Howarth. fc Contil ( 1997 . Figure 2), Figure [5] shows the NIR J — H vs. H — K color- 



color diagram for 609 out of 886 Chandra stars with high-quality (errors < 0.1 mag) JHK photometry listed 



in Ta ble El This plot is valuable for estimating the disk/envelope evolutionary class (I, II, or III, from lLada 



19871 ) of protostars and young stellar objects. Most Chandra sources occupy the color space defined by the 
two leftmost dashed lines (rightward of the locus of main sequence stars, between these two reddening lines 
and including the extension of those lines to the upper right); this region of color space is associated with 
young stars lacking inner disks, referred to as Class III objects, which are reddened by interstellar extinction. 

The ~100 stars to the right of this reddened band are NIR excess sources; most occupy the color space 
between the middle and right-most dashed lines (including sources below the classical T Tauri star locus that 
would fall to the left of the rightmost reddening line if it was extended downward); this space is associated 
with pre-main sequence stars with circumstellar accretion disks, referred to as Class II objects. A handful of 
stars occupy the color space (rightward of the right-most dashed line and its downward extension) associated 
with protostars, referred to as Class I objects. The Herbig AeBe locus shows where these intermediate-mass 
pre-main sequence stars lie in color space; some IR excess sources could be HAeBe stars. The outlier #128 
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is probably an interesting Class I protostar; it is a faint but very hard X-ray source. 

Above the leftmost reddenin g: line is the region of color space occupied by background giant stars 
( Hanson. Howarth, fc Contil 1997 ). Simulations of X-ray detected background giants from the Cepheus B 



study (jGetman et al.1 l2006al ) showed that typically only ~ 3 such objects should be detected in our M17 
Chandra observation; the fact that more than 30 points occupy this region of color space in Figure [5] shows 
that such color-color diagrams are not definitive for classifying sources. NIR photometry can be contaminated 
by diffuse emission or crowding, causing pre-main sequence members of M17 to scatter into this region of 
color space. 

As these classes are based solely on the NIR (JHK) spectral energy distributions (SEDs) of the sources, 
they reflect the evolutionary stage of inner disks. Our analysis is insensitive to the evolutionary status of 
centrally-cleared, outer disks. Thus our Class III status based on JHK colors for the kno wn Q9-B3 stars 



#186, 233, 296, and 488 is not in conflict with the Class I status given to these sources by iNielbock et al 



([20011 ). whose classification is based on IR SEDs that extend out to 20 //m, where outer disks can be detected. 

Figure [6] shows the NIR J vs. J — H color-magnitude diagram for the same stars shown in Figure [5] 
This plot is valuable for estimating individual masses and reddening. Assuming that they lie on the 1 Myr 
isochrone, the majority of ACIS stars appear to be G and F stars (0.5 < M < 2 M Q ) reddened by 3 < Ay < 
15 mag. Sources to the left of the 1 Myr isochrone may be older than the NGC 6618 population ( S|2.5p or 
may have contaminated photometry. 

The drop off of stars towards higher masses is due to the stellar Initial Mass Function (IMF). The 
peak of the mass distribution in the figure is reasonably consistent with the mass completeness limit of 
1.7 M©estimated in ^3.21 for the lightly obscured population. Since the stellar IMF rises steeply towards a 
peak around M ~ 0.3-0.5 M©, most of our X-ray sources are associated with stars having masses below our 
completeness limit. The minimum absorption of Ay ~ 3 mag is consistent with the absorption reported 
in the M17 cluster ( Chini. Elsaesser. fc Neckel 1980 ). Approximately 40 stars are highly obscured with 



Ay > 20 mag; other X-ray discovered members without NIR associations undoubtedly are even more 
obscured. 

Table [6] lists 138 X-ray stars (from the sample of 609 high-quality counterparts) with dereddened mass 
estimates M > 2.0 M©, ordered by mass. Appended to the table are five sources whose photometry has 
large errors but nevertheless strongly suggests that the stars are massive. Some massive stars cataloged in 
the literature are included in this table; they are marked with table notes. Some cataloged massive stars 
(from Tables M and detected by Chandra are missing: ACIS #239, #322, #536, and #543 have no 
photometry available; ACIS #128, #18 2, #399, #647, and #72 6 have dereddened mass estimates less than 
2.0 M©. We adopt the reddening law of Rieke fc Lebofskv ( 1985^ R. Note that inference of stellar masses and 



visual extinction from dereddening involves an ambiguity for stars within the mass range 2-4 M Q (Figure O 
because the reddening vector intercepts the 1 Myr isochrone twice. Table [6] reports the lower of the two 
mass estimates. As a result there is an apparent gap between 4-10 M©in the reported mass distribution, and 
some reported masses in the range of 2-4 M©may be approximately one half their true values. The table 
also reports a NIR (JHK) SED classification for each object (Class I, II, or III described above) obtained 
from the color-color diagram. Of course, these class estimates refer only to hot inner disks and are subject 
to photometric uncertainties, so misclassifications are present. These approximate masses and classifications 



no 



8 The abnormal reddening law of Ry = 4.8 derived for M17 cloud material bv lChini fc Wargaul Jl998h should have little or 
impact on the NIR extinction {Cardelli et al.lll989h : thus we use the standard reddening law from lRieke fc Lebofskv! lll985l) . 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

H-K 



Fig. 5. — NIR color-color diagram. The (light and dark) green circles and red triangles represent sources 
with and without if-band excess, respectively. The dark green circles represent Class I objects, labeled 
by their ACIS sequence numbers. The two blue triangles are foreground stars (#281, #854, ^5.3[) . The 
black solid and gray l ong-dash lines denote the loci of main sequence stars and giants, respectively, from 



Bessell fc Brettl ((1988). The pur ple dash dot t ed lin e is the locus for classical T Tauri stars (CTTS; Class 



II pre-main sequence s tars) from lMever et al.1 (|1997l ). The cyan solid line is the locus for HAeBe stars from 
Lada fc Adamsl ( 1992T ). Three blue dashed lines (marked every 5 mag) are drawn parallel to a reddening 
vector with length Ay = 20 mag, originating from the 05 dwarf, the M4 giant, and the end of the pre-main 
sequence locus. 
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Fig. 6. — NIR J vs. J ~ H color-magnitude diagram. The (light and dark) green circles and red triangles 
represent sources with and without if-band excess, respectively. The dark green circles represent Class I 
objects, labeled by their ACIS sequence numb ers. The two blue tria ngles are foreg r ound stars f H5.3p . The 
purple solid line is the 1 Myr isochrone from ISiess et alJ |2000l) and iBaraffe et all |l998h . The gray dash 
dotted line marks the location of Zero Age Main Sequence (ZAMS) stars; various spectral types are marked 
by black stars. The blue dashed lines represent the standard reddening vectors for various stellar masses, 
marked for every Ay — 5 mag by the blue crosses. Sequence numbers for ACIS sources down to 4 M Q are 
shown in black for cataloged OB stars (from Tabled and in red otherwise. The O star #433 is among those 
omitted due to poor photometry. Sources #51 and #163 are discussed in £)5.3l Typical uncertainties are 
shown graphically; the small median J-band uncertainty (0.02 mag) may be difficult to see. 
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Table 6. Photometrically Selected Intermediate- and High-Mass Stars 
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III 


600 c 


182030 


95- 


161039 


1 


18 


.6 


16.8 


II 


711 


182034 


.85- 


160626 


.2 


19 


6 


16.2 


III 


567 a 


182030 


23- 


161034 


.7 


10 


.6 


16. 1 


III 


1 


181957 


56- 


161204 


3 


9 


.5 


15.3 


III 


650 


182031 


.89- 


161616 


.2 


16 


.6 


14.9 


I 


450 a 


182027 


.80- 


161101 


.4 


12 


.2 


14.6 


II 


201 


182021 


82- 


161123 


.4 


26 


.9 


14.2 


III 


291 


182024 


.33- 


161352 


.8 


17 


.1 


13.6 


III 


347 a 


182025 


.44- 


161115 


.5 


9 


.2 


13.2 


III 


233 a 


182022 


70- 


160833 


6 


6 


.2 


13.2 


III 


413 


182026 


96- 


160347 


.3 


17 


.4 


12.9 


III 


346 


182025 


.42- 


161118 


6 


14 


.7 


12.9 


III 


334 


182025 


08- 


161133 


.7 


18 


,0 


12.9 


II 


875 a 


182101 


07- 


160546 


.7 


3 


.2 


12.8 


III 


296 a 


182024 


.39- 


160843 





7 


.8 


12.7 


III 


118 


182019 


.25- 


161326 


.4 


26 


.9 


11.8 


III 


554 


182030 


03- 


161034 


5 


16 


.2 


11.7 


III 


187 


182021 


.44- 


161206 


6 


19 


4 


11.3 


III 


757 


182037 


.62- 


160332 


3 


21 


.7 


3.7 


III 


27 


182008 


25- 


160552 


5 


7 


.4 


3.7 


III 


691 


182034 


02- 


160800 


6 


23 


.5 


3.7 


II 


817 


182043 


.77- 


161330 


.1 


18 


.7 


3.4 


III 


813 


182043 


.42- 


161245 


9 


5 


.7 


3.4 


III 


488 a 


182028 


66- 


161211 


.3 


3 


.2 


3.4 


III 


466 a 


182028 


.15- 


161049 


.1 


8 


4 


3.4 


III 


608 a 


182031 


.12- 


160929 


6 


7 


.5 


3.3 


III 


582 


182030 


64- 


161028 


3 


16 


.8 


3.3 


III 


560 


182030 


.15- 


161055 


6 


15 


.7 


3.3 


III 


507 a 


182029 


.17- 


160942 


.8 


7 


.5 


3.3 


III 


438 


182027 


62- 


161043 





20 


.9 


3.3 


III 


307 


182024 


56- 


161127 


5 


14 


.7 


3.3 


III 


588 


182030 


.73- 


161101 


.4 


22 


.3 


3.2 


III 


456 


182027 


86- 


160955 


3 


1 


.3 


3.2 


III 


375 


182026 


02- 


161239 


.8 


14 


.8 


3.2 


III 


288 


182024 


.23 — 


161109 


3 


11 


.2 


3.2 


III 


229 


182022 


.58- 


161024 




10 


.9 


3.2 


III 


186 a 


182021 


.43- 


160938 


8 


3 


.9 


3.2 


III 


83 


182016 


.84- 


160729 


.0 


21 


.6 


3.2 


III 


506 


182029 


14- 


161054 


.0 


4 


.1 


3.1 


III 


206 


182021 


92- 


161412 


3 


25 


.5 


3.1 


III 


763 


182038 


.02- 


160310 


.5 


26 


.7 


3.0 


III 


687 


182033 


.76- 


161304 


6 


7 


8 


3.0 


III 


678 


182033 


.21- 


161058 


1 


11 


,0 


3.0 


III 


677 


182033 


17- 


160746 


3 


15 


.4 


3.0 


III 


526 


182029 


.61- 


160603 


.5 


20 


.3 


3.0 


III 


508 


182029 


.20- 


161110 


7 


1 


.6 


3.0 


III 


495 a 


182028 


87- 


161041 


.7 


5 


.8 


3.0 


III 


343 


182025 


.32- 


160939 


.4 


4 


.4 


3.0 


III 


649 


182031 


.86- 


161047 


.1 


20 


,0 


2.9 


III 


638 a 


182031 


70- 


160945 


6 


7 


.0 


2.9 


III 


376 a 


182026 


04- 


161104 


3 


1 


.5 


2.9 


III 


350 a 


182025 


.50- 


161053 


6 


8 


.2 


2.9 


III 


312 a 


182024 


62- 


161108 


4 


6 


.3 


2.9 


III 


253 


182023 


13- 


161131 


.8 


14 


.9 


2.9 


III 


682 


182033 


.43- 


161042 




13 


6 


2.9 


II 


783 


182039 


.71- 


161312 


.5 


7 


8 


2.8 


III 


595 


182030 


89- 


160904 


9 


10 


,0 


2.8 


III 


505 


182029 


.12- 


161247 


.2 


15 


.2 


2.8 


III 


496 a 


182028 


.87- 


161109 


8 


6 


.5 


2.8 


III 


423 


182027 


.20- 


160420 


6 


10 


.8 


2.8 


III 


408 


182026 


82- 


160749 


3 


4 


,0 


2.8 


III 



Table 6 — Continued 



Maws 
("©) 



(est.) 

(5) 



645 


182031 


82- 


161123 


6 


11 


2 


2.7 


564 


182030 


19- 


161213 


4 


17 





2.7 


512 


182029 


32- 


161047 


7 


12 


6 


2.7 


406 a 


182026 


79- 


161057 


8 


7 


4 


2.7 


273 


182023 


76- 


161101 


2 


10 





2.7 


105 


182018 


44- 


161418 


7 


14 


2 


2.7 


324 


182024 


87- 


161127 


2 


6 


2 


2.7 


771 


182038 


82- 


160604 


6 


21 


9 


2.6 


727 


182035 


67- 


160242 


5 


8 


6 


2.6 


607 


182031 


09- 


161125 


9 


17 





2.6 


549 


182029 


95- 


160953 


5 


6 


7 


2.6 


489 


182028 


67- 


160926 





2 


3 


2.6 


179 


182021 


25- 


160944 


1 


15 


8 


2.6 


150 


182020 


32- 


160237 


1 


12 


4 


2.6 


57 


182012 


96- 


161308 


1 


17 


8 


2.6 


800 


182041 


79- 


160241 


7 


9 


5 


2.5 


518 


182029 


44- 


161049 


9 


7 


5 


2.5 


305 


182024 


54- 


160737 


4 


12 


7 


2.5 


710 


182034 


83- 


160750 


9 


9 


3 


2.5 


592 


182030 


80- 


161041 


8 


14 


1 


2.5 


758 


182037 


78- 


160427 


1 


22 


5 


2.4 


568 


182030 


28- 


161131 


2 


8 


6 


2.4 


538 


182029 


82- 


160957 





9 


7 


2.4 


401 


182026 


65- 


160310 


3 


17 


5 


2.4 


30 


182008 


44- 


161409 


4 


8 


2 


2.4 


16 


182004 


72- 


160937 


1 


13 


8 


2.4 


285 


182024 


10- 


160839 


9 


9 


7 


2.4 


628 


182031 


52- 


160910 


1 


9 


6 


2.3 


515 


182029 


39- 


160942 


7 


6 


7 


2.3 


368 


182025 


87- 


160922 


9 


15 


2 


2.3 


236 


182022 


75- 


161123 


2 


9 


7 


2.3 


137 


182019 


95- 


160805 





14 


9 


2.3 


21 


182006 


33- 


160453 


5 


15 


1 


2.3 


537 


182029 


82- 


161057 


2 


12 


7 


2.3 


335 


182025 


10- 


161129 


4 


8 


2 


2.3 


668 


182032 


48- 


161047 





21 


8 


2.2 


655 


182031 


98- 


161030 


2 


14 


1 


2.2 


555 


182030 


04- 


161233 


9 


14 


9 


2.2 


511 


182029 


29- 


161041 


2 


10 


1 


2.2 


437 


182027 


59- 


161133 


3 


10 


3 


2.2 


435 


182027 


48 — 


161013 


2 


3 


4 


2.2 


374 


182026 


01- 


161022 


1 


8 


9 


2.2 


310 


182024 


60- 


161046 


1 


9 


1 


2.2 


264 


182023 


44- 


160953 


1 


9 


9 


2.2 


113 


182019 


06- 


160629 


3 


19 


6 


2.2 


94 


182017 


82- 


160453 


1 


23 


9 


2.2 


680 


182033 


39- 


160918 


6 


6 


3 


2.2 


654 


182031 


97- 


160924 


9 


6 




2.2 


717 


182035 


29- 


161118 


4 


7 


9 


2.1 


694 


182034 


08- 


161044 


3 


4 


6 


2.1 


681 


182033 


41- 


161115 


6 


9 


6 


2.1 


534 


182029 


76- 


160936 


3 


7 





2.1 


434 


182027 


47- 


161358 


4 


4 


7 


2.1 


200 


182021 


79- 


161041 


8 


16 


7 


2.1 


159 


182020 


55- 


160738 





12 


4 


2.1 


142 


182020 


03- 


161330 


8 


17 


9 


2.1 


804 


182042 


04- 


161107 


2 


15 





2.0 


262 


182023 


34- 


160512 


8 


23 


3 


2.0 


225 


182022 


49- 


160951 


5 


10 


7 


2.0 


192 


182021 


62- 


160924 


5 


8 


3 


2.0 


359 


182025 


73- 


161045 


8 


7 





2.0 


7 d 


182000 


65- 


161112 











9 C 


182001 


73- 


160529 











398 b -f 


182026 


62- 


161136 


5 






4: 


433 a -f 


182027 


41- 


161331 









>4 


539 f 


182029 


84- 


161041 


3 






4: 



Note. — Columns 1-2: Source identification from Table fTl 
Columns 3—4: Visual absorption and mass, derived from dcrcd- 
dened location in color-magnitude diagram along the standard in- 
terstellar reddening vector. 

Column 5: Suggested disk/envelope evolutionary class (sec text) 
based on source's position in NIR color-color diagram (Figure 1 5 1 . 
ACI3 #281 and #854 arc likely foreground stars (JOJ and arc 
omitted. ACIS #123 and #543 have inaccurate photometry and 
are omitted. 

a Historically recognized OB stars from Table j*0| 
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b High-mass Class I source from Tablc[s] 

C ACIS #600 is B 163 in Hanson, Howarth, & Conti 11997) and 
is listed in their Tabic 6 as a candidate young stellar object. 

^J-band photometry is not available, however source is bright 
in K-band (10.9 mag) and is likely massive. 

c J-band photometry is not available, however source is bright 
in K-band (9.9 mag) and is likely massive. 

^Photometry has large errors. 
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2.5. X-ray sources not associated with M17 



Approximately 20 extragalactic sources are expected to appear in our ACIS catalog ( £)2.4.1j) . An ad- 
ditional ^20 foreground stars are expe cted to appear i n the c atal og, based on extrap olation of foreground 



contamination simulations presented bv lGetman et alj (|2006a) and I Wang et al.l (|2006f ). This is roughly con- 



sistent with the population of sources observed between the lMy isochrone and the zero-age main sequence 
(ZAMS) track in the color-magnitude diagram discussed in ^2.4.21 (Figure E|. Add itionally, some of those 



sources show large proper motions in the UCAC2 catalog (|Zacharias et al.l 120041 ) and/or are faint, soft, 
minimally-absorbed X-ray sources. Thus, in total ~5% of our catalog of 886 ACIS sources are expected to 
be unrelated to the M17 star forming region. We postpone a detailed membership study until we have a 
longer ACIS observation of M17, where the X-ray spectra will be better characterized and more useful for 
identifying the unrelated populations. 



3. Global Properties of the Stellar Populations 



The COUP study of the Orion Nebula Cluster, whose membership was reliably known in advance 
from many optical and NIR studies, shows a strong statistical association between X -ray luminosity and 
stellar mass in lower mass (M < 3 M Q ) pre- main sequence stars (jPreibisch et al. 2005a). T his relationship is 
confirmed in other young stellar clusters s uch as Cep OB3b and Cep B (jGetman et al.ll2006al) and NGC 2264 



(Flaccomio et al. 2006; Rebull et al 



(Jardine et al 



2006). 



2006 ). The astrophy sical origin of the relation is uncertain: it may arise 



from internal dynam o processes (jPreibisch et al.l l2005af ) or constraints on the external magnetic coronae 



However, the empirical L x — M correlation means that X-ray catalogs should give a surprisingly complete 
census of lightly obscured cloud members d own to a mass limit dependent on the L x limit of the observation 
( Feigelson et al. 20051 Getman et al. 2006a ). Our census of the lightly obscured M17 population is complete 
down to a hard-band X-ray luminosity of logL^ )C ~ 30.5 ergs s " 1 (33. 2|). Assuming that the M17 stars follow 



the same empirical Lh, c — M relation as the Orion population (jPreibisch et al 
that completeness limit roughly corresponds to a mass of 1.7 M 
from non-cluster members ( 32.51) . 



2005a 



Getman et al.ll2006a| ). 



This is achieved with <5% contamination 



A second critical point is that X-ray emission arises primarily from magnetic flaring and thus X- 
ray selected samples are relatively unbiased with respect to the presence or absence of circumstellar disks 
([Feigelson et al.l 120061 ). Some systematic effects are present — Class II stars (systems with accretion disks) 
are systematically ^2 times fainter than Class III stars (systems with weak or absent disks) in the same 
mass s tratum, and a soft X-r ay component from accretion may be occasionally seen — but these are relatively 
small (jPreibisch et al.ll2005ah . Thus the X-ray sample complements IR population studies that often rely on 
the presence of warm dusty disks, although X-ray studies are less complete for the lowest mass objects . From 
the JHK photometric study of the SIRIUS field (somewhat smaller than the ACIS field). I Jiang et al. I (|2002l) 
identified 454 Class I protostars, 2798 Class II /III pre-main sequence stars, and 281 OB candidates. The 
Class I objects are concentrated along the North Bar and South Bar in the Hu region. The Class II/III stars 
are distributed widely across the field with a deficit in the M17-SW region attributable to heavy absorption. 
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3.1. Spatial Distribution of X-ray Stars in M17 

Figure [7] shows the distribution of X-ray stars in the ACIS field, both as individual stars and smoothed 
with a box kernel of width 30". The low obscuration and high obscuration sources are discriminated by their 
spectra as those with median energy below and abov e 2.5 keV respectively . This criterion corresponds to an 



absorption of log Nu ~ 22.3 cm 2 or Ay ~ 10 mag ( Feigelson et al.ll200 



ely. 
3). 



The lightly obscured population shows a simple, centrally concentrated structure centered at (a, 6) = 
(18 ,l 20 m 28 s ,-16°10'50"). This lies about 30" (0.2 pc projected) west of the 04+04 binary CEN 1 (ACIS 
#536 and 543; the field center is shown in §5j). The B2 star CEN 85 (ACIS #466) coincides with the center 
of the cluster distribution. The inner 2' (0.9 pc) radius appears spherical, but the distribution becomes 
elongated towards the northwest around 2'-4' from the center. The central concentration of the cluster 
(surface stellar density > 1 per arcmin 2 ) traced by the ACIS sources is ~8' x 6' (~4 x3 pc). 

Several features are seen in the stellar surface density map for the heavily obscured population (Fig- 
ure Efcd). 



1. The highest concentration of stars (~40) coincides with the main centrally concentrated cluster seen in 
Figure [7h, c. This is undoubtedly the same physical structure exhibiting a wide range of obscuration. 

2. A secondary stellar concentration (~14 stars) is located northwest of the cluster center at 
(18 h 20 TO 22 s , -16°07'30") with 3' (1.4 pc) diameter. This may be part of the northwest elongation seen 
in the lightly obscured population, but there is a distinct dip in stellar density between it and the 
central cluster. This cannot be attributed to absorption: there is a band of CO-emitting molecular 
material west and north of the secondary stellar concentration, but none is seen in the dip region, as 
shown in Figure [TJ). We label this secondary stellar group M17-X in Figure [7]^, where "X" indicates 
X-ray discovered, as it was not previously noted in IR studies. As this concentration is based on only 
14 stars, its reliability is uncertain. 

3. High-density enhancements (~40 stars) in the X-ray star density extend from the main cluster 3' 
towards the west and 4' to the southwest. The densest portion of this enhancement coincides with M17- 
UC1 ( §4.ip and it ends around the KWO ( §4.21) . This arcuate distribution runs along the eastern edge 
of the M 17-SW mole c ular c ore, following a very similar distribution of Class I protostellar candidates 
found by Jiang et al. ( 20021 ). This is an independent detection of this elongated stellar structure; only 



a handful of the ACIS stars are SIRIUS Class I candidates. The ACIS enhancement thus gives a clear 
view of the spatial distribution of the triggered stellar population along the shock front between the 
M17 Hii region and the M17-SW molecular core. 

4. A stellar concentration (10-15 stars) is present at the northern edge of the ACIS field, corresponding 
to the embedded M17-North group of protostars f §4.3|) . 



The spatial distributions of the intermediate- and high-mass stars detected by Chandra are shown in 
Figure [8] based on mass estimates described in §2.4.21 and Table El As expected, the highest concentration 
of massive stars is located around the known bright O stars in NGC 6618. Other concentrations of interest 
include: M17-SW centered on (18' l 20 m 24 s , -16°11'30") around UC1 (§HH), and a tight grouping of three 
stars around the KWO f§4~2|) . 

It is also important to note that several dozen intermediate- and high-mass stars are distributed widely 
throughout the ACIS field. These mostly are lightly obscured and several were identified as possible OB 
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Fig. 7. — Spatial distribution of ACIS X-ray sources. The left-hand panels show the lightly obscured 
population (E mec n an < 2.5 keV or Ay < 10 mag), and the right-hand panels show the heavily obscured 
population (E me dian > 2.5 keV or Ay > 10 mag). The top panels show the individual star positions 
superposed on the Spitzer GLIMPSE 3.6 /im grayscale image, with green contours outlining the North Bar 
and South Bar in the Hn region. The bottom panels show a map of the stellar surface density after smoothing 
with a 30" radius kernel; contours from the MSX 8 /im map are overlaid to show the North Bar and South 
Bar. Four structures discussed in ^3.11 are labeled in panel (d); structure #2 is the newly identified M17-X 
embedded cluster. 
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Fig. 8. — Spatial distribution of X-ray selected intermediate-mass stars (2 < M < 8 M©, small circles) and 
high-mass stars (M > 8 M , diamonds). The MSX 8/xm contours are overlaid. The pluses mark previously 
known O stars ($5] and Table [9]); the known O star not marked as high-mass (missing its diamond) was 
detected (ACIS #433 <^> 01 345) but no mass estimate is reported in Table[6]due to large photometry errors. 
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stars from the early UBV photometric study of lOgura fc Ishidal (|1976l ). Most of these stars do not appear 
to be extremely young. Among the intermediate-mass stars, ~100 have NIR colors consistent with reddened 
Class III objects while about a dozen have if -band excesses characteristic of Herbig Ae/Be stars (Figure[5]). 



3.2. Quantifying the Stellar Population with the XLF 



Establishing and understanding the Initial Mass Function ( IMF) that arises fro m star formation pro- 
cesses are major goals underlying studies of young stellar clusters ( Corbelli et alJl2005h . This effort is subject 
to a variety of difficulties. One of the major problems is the huge level of background and foreground contam- 
ination seen in optical and NIR observations of most clusters. Statistical correction for contami nation using 



observations of nearby control fields or using Galactic stellar models can give uncertain results. iJiang et al 



(|2002f ) used Galactic models of foreground and (extinguished) background stellar populations to correct the 
observed if-band luminosity function (KLF) in regions of the SIRIUS field. The results were satisfactory in 
the limited mass range 0.5 < M < 2 M , but uncertainties in contamination affected the IMF measurement 
at higher masses. For example, Jiang et al. found that it was not possible to derive a corrected KLF for 
region 2C (Figure QJ where the estimated NIR contamination was higher than the number of detected stars. 



X-ray surveys, on the other hand, suffer relatively little contamination ( §2.5j) . and the statistical link 
between X-ray luminosities and masses ($3]) permits a n association between the X- ray Luminosity Function 
(XLF) and the IMF. This i dea was recently pro posed (jFeigelson fe Getman 2005 ) and ap plied to Chandra 
populations in the OMC-1 |Grosso et al.ll2005l ). Cep B/Cep OB3b |Getman et alj|2006ah and NGC 6357 
(jWang et al.l 120061 ) regions. Through their independent analyses, Getman et al. further confirm the tight 
connection between the XLF and the IMF. 

As the total band luminosity systematically underestimates the true emission for obscured objects due 
to the absence of soft band photons, we restrict our analysis to the X-ray luminosities in the hard (2-8) keV 
energy band: the observed luminosity and that corrected for absorption L/j iC . These quantities are 
available for ~600 of the brighter sources in M17 ( §2.3|) . 



Following the approach taken in iGetman et al.l (|2006al ) and in 8j3.ll we divide the 886 Ml 7 ACIS X- 
ray sources into 514 lightly obscured and 372 heavily obscured samples using a threshold of median energy 
Emedian = 2.5 keV, corresponding to Ay ~ 10 mag. Figure [7] shows the spatial distributions of these two 
samples and Figure [9] shows the XLF comparison analysis for these samples using the absorption corrected 
X-ray luminosity L^ jC . For brevity, the similar analysis using Lh is not described in detail. The ~5% 
contamination due to unrelated Galactic or extragalactic sources ( H2.5P is ignored in the following analysis. 

The lightly obscured Orion Nebula Cl uster (ONC) has both t he XLF and IMF acc urate ly measured 



( Fcig elson et al 



2005 



Getman et al.ll2006ah . We assume here, as in IGetman et al.1 (|2006al ) and lWang et al 



2000), that the XLFs of young stellar clusters have the same powerlaw slope at high luminosity and thus the 



ONC can r easonably serve as a ca l ibrato r for more distant clusters. XLFs should be relatively insensitive to 
age effects; iPreibisch fc Feigelsonl (|2005l ) found only a slight (~0.3 in log L t ) decline in luminosity between 
1 and 10 Myr old ONC populations. Specifically, the ONC comparison sample consists of 839 COUP cloud 
members, excluding OB stars and sources with absorbing column density log Nh > 22.0 cm" 2 . The ONC 
sample is nearly com plete down to 0.1 MfT>and about 50% compl ete within the low-mass stellar-sub-stellar 
range of 0.03-0.1 M ? (jPreibisch et alj|2005bl IGetman et al.ll2006ah . 



Luminosities are available for 70% (361 out of 514) and 63% (234 out of 372) of the lightly and heavily 
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Fig. 9. — Hard band X-ray luminosity functions (corrected for absorption) of the lightly obscured (black 
histogram with error bars, left panel) and heavily obscured (right panel) samples of X-ray sources in M17. 
A calibration XLF from the Orion Nebula Cluster is shown as the gray histogram; luminosities below 10 27 
ergs/s occupy the left-most bin. Regression lines to the ONC XLF (gray, solid lines) were scaled upward 
to match the M17 XLFs (gray, dashed lines). The ONC sample is 100% complete down to ~ 0.1 M (black 
dashed line). 



obscured M17 samples, respectively. Uncertainties, AlogL/j , c , are estimated for each luminosity based on 
the number of detecte d counts, using the simulated flux vs. flux uncertainty correlation from Figure 12 of 



Getman et all (2006a) 



Examination of Figure [9] shows that the ONC and M17 XLFs are nearly parallel at the high luminosity 
end, and the M17 source counts turn over due to completeness limits around logLh.c ~ 30.5 ergs s _1 
(M ~ 1.7 Mq) for the lightly obscured sample and log£h., c ~30.9 ergs s" 1 (M ~ 2.6 M Q ) for the heavily 
obscured sample. We thus scale the ONC XLF regression line vertically upward to estimate the total 
population of M17 in the ACIS field. Due to systematic uncertainties, the XLF-L^ and XLF-L/ l c analyses 
provide slightly different estimates for the total population of the lightly obscured sample in M17 down to 
0.1 Mq, so we report a range for this quantity of 2500 — 3500 stars. About 60% of this population lies in the 
concentrated NGC 6618 cluster. Similar uncertainty for the heavily obscured population gives an estimated 
range of 5000 — 7000 obscured objects. Thus the total estimated population of young objects in M17 down 
to 0.1 Mp, within the 17' x 17' ACIS-I field is 7500-10500 stars. 
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4. X-rays from embedded populations 



4.1. The M17-UC1 region 

Radio and IR studies have shown that star formation around the M17 Hn region is most active in the 
M17 -SW cloud. The best-studied subregion is around the massive binary IRS 5N and 5S with separation 5 " 



(e.g. iFelli et alJEiii Ijohnson. Depree. fc Gos]ll998l : Ichini. Nielbock. fc Becklboool : iNielbock et alll200lh . 



IRS 5N has bolometric luminosity Li> i ~ 5000 L© peaking in the MIR and ionizes the rapidly evolving 
ultracompact Hii region M17-UC1 <^> G15.04-0.68. IRS 5S with L hoi ~ 2000 L© is radio-quiet and less 
absorbed. They appear to be massive Class I protostars, probably early-B stars. The interstellar matter 
around these central stars is complex with dense dusty cloud clumps and an arc-shaped ionization front. 

The top panels of Figure [TO] show the deep K and Chandra images of IRS 5 and its vicinity. ACIS #322 
did not arise from the source finding process but was manually added to the catalog at the known position 
of the MIR source IRS 5N <^ UC 1. Although its formal detection significance (column 13 of Table [2|) is low 
(0.9), confidence that X-ray emission is seen at that position is bolstered by the fact that the 5 nearby X-ray 
photons are harder (E me< a an = 4.7 keV) than the local background (Emedian — 3.5 keV). IRS 5S (#309) 
is definitely present with 19 photons. Its X-ray spectrum indicates absorption of log Nh ~ 22.8 cm -2 and 

Thi s source was recently 
20051 k the X-ray source is 



intrinsic (absorption-corrected) hard band luminosity logL/j iC ~ 3 0.9 ergs s 



resolved into 3 components using a high-resolution if-band image (jChini et al 
spatially associated with the brightest X-band component. 

A cluster of ^20 X-ray stars is present within 20" (0.15 pc) radius of the IRS 5 binary, most with NIR 
counterparts. These sources are heavily absorbed and have intrinsic X-ray luminosities near the top of the 
XLF. Source #324 is associated wit h the NIR star B 279 (K = 10.75 mag) which is likely responsible for 
powerful methanol and OH masers (jNorris et al.l I1987H . It emits logL^ iC ~ 30.9 ergs s _1 with absorption 
\ogN[j ~ 22.4 cm -2 . Source #334 is associate d with the IR sta r B 273 with K = 10.3 mag and strong 
CO bandhead absorption but no emission lines ( Chini et al. 20051) . It emits \ogLh, c ~ 30.5 ergs s _1 with 
log Nh ~ 22.4 cm~ 2 . Chini et al. interpret this source as an FU Orionis type object, a massive young star 
with heavy accretion. 

This concentration of ~20 X-ray luminous stars around IRS 5 has about 1/5 the number of similarly 
luminous stars seen in the lightly obscured NGC 6618 cluster fi j3.2[) . implying that the total population of 
the embedded cluster is ^300 stars. Source #329, one of the weaker X-ray sources in the region, showed 
flaring behavior and is probably one of these lower mass stars. However, this inference is quite uncertain, 
and inconsistent with the relatively small population found in the Becklin-Neugebauer region behind the 
Orion Nebula which is similar to the M17-SW IRS 5/UC1 complex. In Orion, the deep COUP image shows 



only a sparse cluster of lower mass stars accompanying the IR-luminous group of proto-B stars (jGrosso et al 



2005). We look forward to more directly measuring the population around IRS 5 in an upcoming deeper 
Chandra observation. 



A group of dense molecular cores ( Peak 3 in the survey of Vallee fc Bastienl Il996l ) and masers are 
concentrated about 30" west of M17-UC1. Johnson. Depree. fc Goss ( 19981) argue that this region is the one 
most recently triggered by the shock from the main M17 Hn region. Two ACIS source s (#239 and #244) are 
positionally coincident with the luminous Class I IR source Anon 1 ( §4.4)1 . reported bv lNielbock et al.l (|200lh 
to have an extremely steep spectral index. The properties of source #239 suggest it is the more likely true 
association; it is unusually X-ray luminous, constant, and heavily absorbed (logLh jC ~ 31.7 ergs s _1 and 
log Nh ~ 23.5 cm -2 ), and has no identified NIR or MIR counterpart. The observed absorption is roughly 
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Fig. 10. — M17-UC1 (upper panels, ^4. 1|) and the Kleinmann- Wright Object (lower panels, £14. 2p seen in the 
SIRIUS if -band (left) and Chandra X-ray (right) images. The X-ray gray scales span 1 to 21 counts (upper 
right) and 1 to 9 counts (lower right) in 0.5" x 0.5" pixels. ACIS sources are marked by diamonds (left) 
and sequence numbers (right); source extraction regions matching the Chandra Point Spread Function are 
shown in blue. SIRIUS counterparts (+) range from 9 to 14 mag (upper left) and from 7.5 to 15.5 mag 
(lower left) in K. The radio arc of ionized gas associated with M17-UC1 can been faintly seen in the if-band 
image (upper left). ACIS #309 matches the brightest IR component of IRS 5S. 
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half of the estimated total absorption through Peak 3 ( Hobson et al. 1993h . consistent with the star lying 
at the center of the core. The second ACIS source, #244, shows flare-like variability (Figure [3]) and more 
typical properties (logLft )C ~ 31.1 ergs s _1 , log A# ~ 22.2 cm -2 , K ~ 12 mag counterpart), suggestive of a 
young star lying in front of the star forming cores though still inside the M17 giant molecular cloud. 



4.2. The Kleinman- Wright Object and its cluster 



The lower panels of Figure [10] show the deep K and C handra images around the IR -luminous star known 
as the Kleinmann- Wright Object (KWO) or M17-IRS1 (jKleinmann fc Wrightl Il973l) . IR study reveals a 
photosphere with T e ff = 25,000 K and L^ i — 5 x 10 3 L© with strong variable hydrogen emission lines 
(jChini et al.ll2004aT ). It is modeled as a Hcrbig Be star of spectral type BOV lying inside a dusty envelope 



with M ~ 10 M Q , r ~ 0.1 pc, and Nh ~ 2 x 10 cm . It may be the most heavily absorbed and 
youngest Herbig Be star known, with a Class I protostellar IR spectral energy distribution. Chini et al. 
resolve a nearby bright NIR source (KWO component 2) lying 1.2" northwest of the Herbig Be star (KWO 
component 1). Chini et al. further find a concentrated cluster of about 150 heavily reddened JHK stars 
within 0.7 pc (70") of the KWO. 

ACIS source #123, with 31 counts lies 1.2" n o rthwes t of SIRIUS source 15128; SIRIUS does not 
resolve the two components described bv IChini et alJ (|2004al ). although the SIRIUS J-band image appears 
elongated to the northwest, in the direction of the ACIS source. Comparing the SIRIUS and ACIS positions 
to components 1 and 2 of the KWO, based on positions obtained by Chini et al. in the z-band with the 
ESO NTT telescope, it appears that the SIRIUS source closely matches KWO component 1 (the Herbig Be 
star, which is brighter in the NIR), while ACIS #123 closely matches KWO component 2. The Chini et 
al. positions are systematically offset from the ACIS position and the SIRIUS position (after registration 
to ACIS/2M ASS) by ~ 1 " to th e southwest. The KWO Herbig Be star is not detected in this Chandra 



observation. IChini et al.l (|2004al ) report that component 2 (o- ACIS #123) is probably a member of the 
KWO cluster. It exhibits \ogLh, c ~ 31.0 ergs s _1 and log Nh ~ 22.5 cm -2 equivalent to Ay ~ 20 mag. 
The situation resembles that of the Becklin-Neugebauer Object in th e OMC-1 where a p reviously unknown 
lower-mass companion greatly outshines the young B star in X-rays ( Grosso et al.ll2005l ). 



A dozen other ACIS sources are present within 30" of the KWO, a higher density than generally seen 
in the M17-SW cloud though less rich than seen around IRS 5/UC 1 ( ij4.ip . Most of these sources have 
11 < K < 15 mag counterparts consistent with intermediate-mass stars in the embedded cluster. As with 
UC1, an upcoming deeper Chandra observation should detect more members of the cluster around the KWO 
and perhaps the Herbig Be star itself. 



4.3. The M17-North cloud 



The M17-North cloud cores are dusty molecular con centrations extending over 2' at a location around 10' 
north of the bright M17 Hn region ( Wilson et al. 2003). Star formation here is probably occurring without 
triggering by the Hn region. It has five em bedded IR source s; radio and ionized line emission suggest that 
IRS 3 is a late-0 star behind the cloud core ((Klein et alJll999h . IRS 1 and 3 are off the ACIS field; IRS 2 is at 
the field edge; IRS 4 and 5 are not detected in the Chandra image. However, the ACIS image (Figure [2]) and 
the smoothed map of obscured ACIS sources (Figure right panel) show a density enhancement of a dozen 
absorbed sources. For example, source #618 lying 10" from M17-North IRS 4 has \ogLh ~ 30.8 ergs s" 1 
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and log Njj ~ 22.5 cm 2 without any known counterpart. Two sources (#660 and 665) lie 6" from each 
other, again without if-band counterparts. 

Assuming a standard XLF ( i)3.21 Figure O, the presence of this group of X-ray luminous stars implies 
that dozens or hundreds of lower mass stars are present. As with the IRS 5/UC1 and KWO clusters, the 
population of the clusters cannot be reliably discerned from the present shallow X-ray exposure. However, 
we note that the (probably untriggered) M17-North group appears to be spread over a larger region (~2'or 
1 pc) than the (probably triggered) clusters in the M17-SW cloud. 



4.4. Protostar candidates 

Chandra studies of nearby low-mass star forming regions (such as the Ophiuchus, Taurus, and Serpens 
clouds) have established that Class I protos tars are X-ray emitters with typical luminos i ties of ' 



flaring light curves, and hard spectra (e.g. Ilmanishi et al.l l2001u iGetman et al 



case for X-ray emission from Class sources is problematic (jTsuiimoto et al 



TO 30 ergs s 1 , 



2005a: Prcibisch 2004). The 



2005; Getman et al 



2006bl ). 



Systematic studies of X-ray protostars in MSFRs, however, have not yet emerged. 

From the JHK study of Jiang et al. ( 20021 ). 157 candidate Class I protostars were reliably identified 
from strong if-band photometric excesses. About half lie in the North Bar and South Bar in the Hn region, 
while the other half are widely distributed across their field. Others can be identified by MIR excesses in 
the GLIMPSE survey. 

Table [7] lists the 64 Chandra sources associated with protostellar candidates satisfying one or more of 
the following criteria: 



1. Strong if-band excess in the (J — H) vs. (H — K) diagram beyond the locus of reddened Class II T 
Tauri stars (Figure [5]). 

2. Excess in the shortest wavelength band of the Spitzer IRAC detector. Observations of Taurus-Auriga 
young stars with reliable classifications show that the color criterion K — [3.61 > 1 -5 mag isolates Class 
0-1 protostellar systems from Class II-III T Tauri systems (jHartmann et al.l 120051 ). 



3. Exces s in the second IRAC band with [3.6] — [4.5] > 0.7 mag also separates protostars from T Tauri 
stars ( Hartmann et al. 20051 ). 



Uncertainties on the K — [3.6] and [3.6] — [4.5] colors were calculated, and only sources exceeding the stated 
color thresholds by more than 1 standard deviation are reported. X-ray luminosities in the hard band are 
reproduced from Table [3] when available; for the fainter sources, we approximate logL/j ~ 29 ergs s _1 , based 
on the ratio between counts and Lh for heavily absorbed brighter stars in Table [3] The luminosities range 
up to 10 310 ergs s -1 , similar to the distribution of luminosities of Class I sources in the p Ophiuchi cloud 



(jlmanishi et al.ll200ll ) 



Only 13 out of 157 Class I candidates identified by i^-band excess (criterion 1 above) are detected. The 
paucity of extreme if-band excess sources can also be seen in the NIR color-color diagram discussed in £12.4.21 
(Figure[5]). The median K magnitude (~13.5) for the 13 X-ray detected sources is 1 magnitude brighter than 
the median for the 157 Class I NIR candidates, suggesting that, at the sensitivity limits of our observation, 
only the higher-mass end of the protostar XLF is detected. An upcoming longer Chandra observation should 
detect lower-mass Class I candidates. 
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An additional 51 Class I candidates in Table [7] identified by colors from longer wavebands (criteria 2 and 
3 above), more sensitive to circumstellar disks, are detected. Many of these are highly obscured; half have 
Emedian > 3 keV (corresponding to Ay > 15 mag). These 51 sources (and perhaps some of the ~80 X-ray 
sources without counterparts described in ij2.4.ip represent a significant population of Class I candidates 
that were not previously available from NIR studies. 
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Table 7. Photometrically selected X-ray emitting candidate protostars 



Seq 



CXOU 



K 



log N H 



E„ 



Photometric Criteria 









(mag) 


(erg s x ) 


(cm- 2 ) 


(keV 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




182000.65- 


161112.0 


10.88 


30.40 


22.0 


2.0 


21 a 


182006.33- 


160453.7 


12.52 


29: 




3.3 


26 


182008.21- 


161040.9 


13.72 


30.33 


22.9 


2.8 


30 a 


182008.44- 


161409.6 


11.62 


30.71 


21.8 


1.9 


47 


182011.77- 


160512.2 


14.35 


29: 




2.3 


51 a 


182012.31- 


160447.7 


10.35 


30.27 


22.3 


2.3 


53 


182012.44- 


160610.0 


15.73 


29: 




1.1 


57 a 


182012.96- 


161308.4 


12.46 


29: 




2.0 


58 


182013.09- 


161238.1 




29: 




3.4 


60 


182013.48- 


160912.2 


14.45 


29: 




1.8 


74 


182015.80- 


161033.2 


14.99 


29: 




2.3 


81 


182016.54- 


161003.0 


13.83 


30.99 


23.1 


4.1 


94 a 


182017.82- 


160453.4 


13.67 


31.11 


22.3 


3.3 


95 


182017.82- 


161401.2 


13.75 


30.79 


22.8 


3.7 


105 a 


182018.44- 


161419.0 


12.22 


30.52 


22.6 


2.0 


106 


182018.51- 


161421.4 


12.76 


30.38 


22.2 


2.6 


108 


182018.65- 


160607.8 


14.07 


29: 




4.1 


109 


182018.83- 


161424.8 




30.45 


22.8 


3.6 


113 a 


182019.06- 


160629.6 


12.96 


30.84 


22.9 


3.1 


117 


182019.24- 


160800.3 


13.01 


30.44 


22.5 


2.7 


128 


182019.60- 


161038.8 


11.92 


30.45 


23.2 


4.8 


137 a 


182019.95- 


160805.2 


12.44 


29: 




3.1 


156 


182020.53- 


160456.3 


14.15 


30.84 


22.6 


3.1 


179 a 


182021.25- 


160944.3 


12.45 


29: 




2.5 


200 a 


182021.79- 


161042.0 


12.81 


30.72 


22.6 


2.6 


201 a 


182021.82- 


161123.6 


11.36 


29: 




4.4 


207 


182021.97- 


161100.6 


13.31 


29: 




3.1 


209 


182022.01- 


161024.8 


14.06 


29.99 


22.3 


2.0 


221 


182022.40- 


161403.9 


13.29 


29: 




3.2 


234 


182022.70- 


160951.5 


13.88 


30.36 


22.6 


2.7 


235 


182022.70- 


161613.8 


14.90 


29: 




3.5 


304 


182024.51- 


160644.3 


15.06 


30.56 


22.7 


2.7 


305 a 


182024.54- 


160737.6 


12.00 


30.18 


22.2 


3.7 


324 a& 


182024.87- 


161127.5 


10.75 


30.74 


22.4 


2.7 


328 


182024.93- 


160459.5 


14.38 


29: 




3.2 


363 b 


182025.82- 


161121.0 


13.61 


30.32 


22.2 


2.1 


372 


182025.95- 


160938.2 


13.63 


29: 




2.5 


380 


182026.12- 


160455.9 


14.24 


31.07 


22.8 


3.6 


438 a 


182027.62- 


161043.3 


12.25 


29: 




3.8 


472 


182028.28- 


161130.5 


11.08 


30.52 


22.4 


2.6 



JH,HK 

(7) 



K-[3.6] 
(8) 



[3.6]-[4.5] 
(9) 



V 
V 



V 
V 



V 
V 
V 
V 
V 
V 

V 

V 
V 
V 
V 
V 
V 
V 
V 

V 
V 

V 
V 
V 
V 
V 



V 

V 
V 
V 

V 

V 
V 
V 
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The two sources from Table [6] identified as Class I (#324 and #650) are two of the if -band excess 
sources in Table [7J An additional 21 sources in Table [7] also appear in Table [H Assuming that the protostar 
criteria defined above supersede the estimated disk/envelope evolutionary classes in Table [51 these 23 sources 
(marked by note a in TableUJ define a new sample of X-ray selected intermediate- and high-mass protostars; 
none are known OB stars. These stars especially warrant further study. 



Nielbock et al.l (|20011 ) report an independent survey for high-mass Class I stars in M17 based on ground- 
based 10 fim. and 20 fim imaging. Of their 22 luminous IR sources, we detect the 9 listed in Table [5] (nominally 



Table 7— Continued 





PXOTT 


K 


log Lh 


log iVf/ 


K 

^median 


Photometric Criteria 








\ erg b ) 


(cm ) 


^ivc V J 


JH,HK K-[3.6] 


[3.6]-[4.5] 


UJ 






(A\ 


\°) 




(7) (8) 


(9) 


498 


182028.95-160957.7 


13.26 


30.43 


22.4 


2.2 


V 




500 


182029.01-160831.5 


14.52 


29: 




2.9 


V 




535 


182029.80-161125.4 


12.46 


29.72 


22.2 


1.8 


/ 

V 




562 


182030.18-161216.5 


14.15 


29.86 


22.4 


2.3 


V 




595 a 


182030.88-160905.1 


11.47 


30.55 


21.9 


2.1 


V 




596 


182030.90-160954.7 


12.69 


29: 




2.3 


V 




607 a 


182031.09-161126.1 


12.42 


30.43 


22.0 


2.6 


V 




610 


182031.23-160218.2 


12.71 


30.83 


22.7 


3.1 


V 




637 


182031.70-160612.1 


14.64 


29: 




4.3 


V 




643 


182031.79-160349.3 


14.26 


29: 




2.9 


V 




650 a 


182031.89-161616.5 


9.57 


30.07 


20.6 


1.7 


V V 




652 


182031.94-160239.7 


13.13 


31.30 


22.8 


3.8 


V 




654 a 


182031.97-160925.1 


11.26 


30.50 


21.9 


2.6 


V 




668 a 


182032.47-161047.2 


13.54 


29: 




4.4 


V 




677 a 


182033.17-160746.6 


12.10 


29: 




2.7 


V 




679 


182033.32-161050.9 


13.63 


29: 




3.1 


V 




682 a 


182033.42-161042.7 


11.79 


29: 




3.4 


V 




709 


182034.79-161053.8 


11.66 


30.47 


21.8 


1.9 


V 




718 


182035.34-161325.4 


13.58 


29: 




2.3 


V 




757 a 


182037.62-160332.6 


12.00 


30.95 


22.3 


2.7 


V 


V 


762 


182037.99-160925.8 


12.78 


30.62 


21.8 


2.0 


V 




789 


182040.93-161111.7 


15.67 


29: 




1.3 


V 




782 


182039.62-161146.5 


12.41 


30.70 


22.8 


3.1 


V 




784 


182039.83-160500.1 


15.41 


30.46 


22.8 


3.3 


V 





Not e. — Columns 1—2: Source identification from Table [T] Column 3: /("-band magnitudes from I Jiang et al 
(2002). Columns 4—5: X-ray spectroscopic parameters from Table [3] when available. For faint sources, a value 
of logih ~ 29 ergs s _1 is assumed (see text). Column 6: Median X-ray energy from Table [T] Columns 7—9: 
Selection criteria from H4.4I 

a Source is also classified as an intermediate- or high-mass star (Table |6}. 

b These protostars lie in the M17-SW IRS 5 cluster around the ultracompact Hll region UC1. 
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including the tentative Chandra detection of IRS 5N) . The X-ray non-detections are not systematically fainter 
in the IR or more obscured. Four of the sources in Tableware cataloged OB stars (ACIS #186, 233, 296, and 
488), but the remaining five stars (ACIS #128, 239, 309, 322, and 398) should be added to the 23 described 
above as X-ray detected intermediate- and high-mass protostar candidates. Only the first of those, #128, 
also appears in our list of candidate protostars (Table [7]). However it is not identified (in Table [6]) as massive 
based on its NIR colors; it is a distinct outlier in the JHK color-color diagram (Figure [5]) and is a hard 
X-ray source. Source #239, has no NIR counterpart; #309 and #398 appear in Table [6] but not in Table [7] 
Source #322(<=> IRS 5S; JO}, has a SIRIUS counterpart but unreliable JHK photometry. 

All the sources in Table [8] have X-ray luminosities arou nd 10 30 ergs s" 1 . T his is typical of late-0 and 
early-B stars without strong radiatively accelerated winds (jStelzer et al.l 120051 ) . In some cases, the X-ray 
emission probably arises from unresolved low mass companions. None of these 9 stars are powerful wind 
sources like the central 04+04 binary ( §5.2| . A useful comparison can be made to the cluster of proto-B 
stars in the Becklin-Neugebauer region of the Orion Molecular Cloud 1 (OMC-1). In OM C-1, one star 
(Sour ce n) has logL^ > 30 ergs s _1 and three others are detected around logL^ ~ 29 ergs s _1 (jGrosso et al. 
2005). In M17, there are a dozen stars similar to Source n. These stars are not confined to the well-studied 



M17-UC1 and KWO clusters; indeed, several appear on the western side of the image away from the M17 
Hn region. This supports the claim that star formation in the M17 cloud is w idely distributed, as indicated 
by the widespread distribution of JHK protostars seen by I Jiang et all |2002h . 



Finally, three interesting high-m ass sources merit comment. Source #488 IRS 15) has recently 
been described bv lChini et al.l ([20061 ) as a newly- formed 26 M star (spectral type B0.5V) that has stopped 
accreting but is still surrounded by a huge remnant disk. Its X-ray spectral properties are similar to other 
early-B stars in M17. No X-ray source is found associated with M17-S01 (M17 silho uette object 1), an 
intermediate- or high-mass star with a dusty envelope and bipolar reflection nebula ( Chini et al. 2004bl 
Sako et al.ll2005af ). CEN 92, a r eddened, massive binary system with Lm > 10 5 and a circumstellar 



IR-bright shell ( Chini et al1l2005l ). is also undetected in X-rays. Chini et al. note that CEN 92 and IRS 5N 



(ionizing M17-UC1) might be similar objects; both are emission-line stars and IRS 5N might be binary. 



4.5. Distributed star formation across the molecular cloud 



A number of heavily obscured X-ray stars lie in other dense molecular regions of M17-SW. IVallee fc Bastien 



(|l996l ) identify six peaks in dust and molecular emission maps as dense cores with sizes 10 — 30" (0.1 — 0.3 pc), 
masses 70 — 470 M Q , densities 1 — 8 x 10 5 cm~ 3 and absorbing columns of 1.5 — 3 x 10 23 cm~ 2 . ACIS sources 
#236 in Peak 1, #97 in Peak 2, #177 in Peak 4, and #275 in Peak 5 all have intrinsic X-ray luminosities 
around logLh iC ~ 30.5 — 31 ergs s _1 and absorptions around log Njj ~ 22.3 — 22.6 cm" 2 . The luminous 
source #239 44> Anon 1 in Peak 3 was discussed earlier in §4.11 and Table El None of these molecular gas 
concentrations appears to have produced rich stellar clusters, but the high X-ray luminosities of the detected 
stars implies that more sources with lower luminosities are present. 

We can examine the X-ray source list for heavily obscured stars that have not been studied at other 
wavelengths. Due to the short exposure, only the stars with logL/j > 30.5 ergs s" 1 are sufficiently bright 
to permit quantitative estimation of high levels of obscuration. Six ACIS sources have spectral fits with 
log N H > 23.3 cnT 2 equivalent to Ay > 100 mag: #37, 135, 239, 246, 255, and 706. Two of these (#239 
and 246) reside in the M17-UC1 cluster while the others are not in known molecular cores. 



Altogether, we support the evidence of I Jiang et al.l (|2002h and others that star formation is widely 



Table 8. X-ray emitting high-mass Class I sources from 



Nielbock et al, 



(|200lh 



Seq 


cxou 


Name 


K 


OLK,N 




SpTy 


logL h 


log7V ff 








(mag) 




(L ) 




(erg s -1 ) 


(cm -2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


128 


182019.60-161038.8 


B 353 


11.92 


2.1 


75 




30.45 


23.2 


186 


182021.43-160939.1 


CEN 35 


10.42 


0.8 


<140 


B3 


30.42 


22.1 


233 


182022.69-160833.9 


CEN 16 a 


9.10 


1.3 


<275 


B0: 


30.69 


21.7 


239 


182022.87-161148.5 


Anon 1 




3.3 


590 




30.78 


23.5 


296 


182024.39-160843.3 


CEN 31 a 


9.40 


2.3 


910 


09.5 


29.58 


21.6 


309 


182024.60-161139.2 


IRS 5S 


9.34 


2.5 


1900 




30.50 


22.8 


322 


182024.83-161135.3 


IRS 5N b 




4.2 


4775 




29: 




398 


182026.62-161136.5 


IRS 10 


10.56 


1.5 


55 




29.99 


22.5 


488 


182028.65-161211.6 


IRS 15 a 


10.08 


3.0 


1405 


B0.5 


30.78 


22.2 



Note. — Co l umns 1—2: Source identification from Table [TJ Column 3: Counte rpart identification from 
Nielbock et aL ( 2001). Column 4: if-band magnitudes from I Jiang et ah J 2002). Column 5: Spectral 



index from iNielbock et alJ ( 200ll ). Column 6: Infrared luminosity from iNielbock et al.l ( 200ll ). Column 
7: Spectral type. Columns 8—9: X-ray spectroscopic parameters from Table [3J 

a These stars are surrounded by resolved infrared-brig ht dusty disks |Chini et alj|2005l bood ). 

b Sourcc =#=322 is tentative (i M.ip and its luminosity is assumed. 
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distributed across the cloud in addition to concentrations around M17-UC1, the KWO, and M17-North. 
This implies either that star formation triggered by the passage of the ionization front of the large M17 Hu 
region does not always produce concentrated clusters and/or that much of the star formation in the cloud 
was not triggered. 



5. X-rays from massive stars 



Individual O stars have long been known to emit soft X 



(Chlebowski et al. 1989; Berghoefer et al. 1997; Sana et al 



r ays a t levels L x oc 10 7 Lboi ~ 10 31 33 ergs s 1 
2006) , arising from a myriad of weak shocks 



withi n their radiatively accelerated winds (|Lucv fc White! ll98*oT iPallavicini et al.l 1 1 9 8 it lOwocki fe Cohen 
1999). While this soft (kT ~ 0.5 keV), slowly varying X-ray component dominated early studies of 



O stars, recent studies reve al discrepancies with X-ray line widths, inferred densities, and temperatures 
(jWaldron fe Cassinelrill200ll ). In particular, a moderately hard componen t (kT ~ 2 — 3 keV) may be present 
which is sometimes rotationally modulated and / or shows rapid variability ( Feigelson et al .1120021 : Schulz et al 



2003tlstelzer et alj|2005t Ischulz et allhood) . While not fully understood, this component has been attributed 
to a large-scale magnetically confined (or at least c hanneled) wind shock that diverts radia l outfl ow into an 
equatorial disk, causing an X-ray emitting shock (jBabel fc Montmerld Il997[ iGagne et al.l I2005T ) . A lumi- 



nous hard X-ray component, sometimes reaching Ly t ~ 10 ergs s and kT > 10 keV, can be produced 
by wind-wind collisions in close massive binaries, particula rly if one component is a Wolf-Rayet star (e.g 



Skinner et al 



2002 



Pollock et al. 2005 : Skinner et al. 20061) . or possibly by a non-thermal inverse Compton 



process (IChen fc White 1991). An emerging class of l ate-O/early-B emission line stars known as "7 Cas 
analogs" ( Smith fc Balona 20061 Rakowski et al. 2006 ) show very hard X-ray spectra (kT > 8 keV) and 
variable X-ray ligh tcurves, perhaps due to some combination of winds, magnetic fields, and disk processes 



(jSmith et al.l 12004). Wind-generated X-rays become insignificant in stars cooler than early-B types, but 
these systems are sometimes detec ted in X-rays due to lower mass pre-main sequen ce companion s (e 



Stelzer et al.ll2005t ISana et al.ll2006l ) or perhaps due to their own intrinsic emission (e.g. IStelzer et al]|200' 



M17 provides an excellent laboratory for study of these OB wind-generated X-rays due to the large 
population of coeval and codistant massive stars in a single field. The main limitation is that we cannot 
uniformly measure their soft X-ray components due to heavy absorption towards some systems. 

Table [9] gives t he X-ray properti es of all published O and early B stars in the region (from TFM03, with 
IRS 15 added from IChini etalll2006h Fl. The stars are listed in order of decreasing mass; see the table notes 
for descriptions of the tabulated quantities. While all of the O stars in Tableware detected with Chandra, 
only 19 of the 34 B0-B3 stars are detected. 

Figure [Tl] shows a close-up view of the cluster center where many of the OB stars are concentrated. 
Figure [12] shows the X-ray spectra of the components of the 04+04 binary and four other OB stars, 
illustrating the diversity of plasma properties. Note that all but one of these spectra show very hard thermal 
plasma components with kT > 4 keV. As mentioned above, s uch hard X-ray emiss ion is most commonly 
seen in colliding wind binaries; a famous example is rj Carinae ([Corcoran et alj |2004). 



9 The two O stars listed in TFM03 as Chandra non-detections (CEN 102 and CEN 34) were determined by 
lHanson, Howarth, fc Contil lll997h to be background red giants based on their NIR brightness and the presence of CO bandhead 
absorption in their IR spectra. The B stars CEN 95 and CEN 33, also not detected in X-rays, were similarly identified by 
Hanson et al. to be background sources. These four sources are omitted from Table [9] 
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Table 9. X-ray properties of cataloged OB stars in NGC 6618 



Optical Properties 



X-ray Properties 



(mag) 
(4) 



< L o> 

(5) 



# 
(6) 



(") 
(7) 



kT 
(kcV) 
(10) 



(erg s - 
(11) 



(erg a" 1 ) 
(12) 



CEN 43 


03-04 


10396 


8.18 


5.8 


675 


0.2 


82 


22.1 


2.0 


31.02 


31.53 


CEN la 


04 




6.59: 


5.7 


536 




1912 


22.4 


0.9 + 15 


32.37 


33.16 


CEN lb 


04 




6.59: 


5.7 


543 




3866 


22.3 


0.7+10.4 


32.78 


33.28 


CEN 2 


OS 


10399 


7.49 


5.5 


701 


0.1 


250 


21.2 


1.6 


31.14 


31.69 


CEN 37 


03-06 


10398 


7.88 


5.5 


574 


0.2 


90 


22.4 


0.6 


30.42 


32.25 


OI 345 


06 


14703 


9.18 


5.3 


433 


0.0 


263 


22.0 


0.6 


30.32 


32.06 


CEN 18 


07-08 


28446 


7.84 


5.1 


366 


0.2 


12S 


21.8 


4 . 4 


31.02 


31.31 


CEN 25 


07-08 


7877 


8.94 


5.1 


593 


0.1 


16 


21.9 


0.7 


29.53 


31.00 


OI 352 


OS 


... b 


6.82 


5.0 


732 




1375 


22.2 


0.5 + 4.4 


32.00 


32.95 


CEN 16 


09-B2 


28651 


9.10 


4.8 


233 


0.4 


73 


21.7 


3.5 


30.69 


31.01 


CEN 61 


09-B2 


7833 


9.23 


4.8 


567 


0.1 


44 


22.3 


0.6 


30.05 


31.72 


CEN 3 


OO 


10401 


7.59 


4.8 


720 


0.2 


108 


21.6 


0.6 


29.72 


31.34 


OI 174 


09 




7.80 


4.8 


12 


0.5 


121 


22.1 


0.5 


30.10 


32.04 


CEN 31 


09.5 


27457 


9.40 


4.7 


296 


0.1 


10 


22.0 


2: 


28.81 


30.31 


CEN 92 


B0 


14714 


9.10 


4.5 












<29.0 


<29.5 


CEN 28 


B0 


8849 


11.40 


4.5 


647 


0.1 


4 


21: 


2: 


28.5: 


28.8: 


OI 582 


B0 


... d 




4.5 


875 




19 


20: 


1.7 


29.69 


30.20 


IRS 15 


B0. 5 


5988 


10.08 


4.3 


488 


0.0 


71 


22.2 


2.2 


30.78 


31.30 


CEN 57 


Bl 


14554 


9.58 


4.0 


347 


0.2 


117 


22.1 


2.5 


31.00 


31.44 


CEN 101 


Bl 


9090 


11.39 


4.0 


496 


0.0 


111 


22.0 


1.8 


30.76 


31.32 


CEN 97 


Bl 


13963 


10.78 


4.0 


507 


0.0 


134 


22.2 


2.9 


31.11 


31.52 


CEN 100 


Bl 


17147 


10.85 


4.0 


608 


0.0 


223 


22.1 


3.9 


31.41 


31.75 


CEN 45 


Bl 


4442 


10.81 


4.0 


726 


0.2 


11 


21.8 


2: 


29.71 


30.19 


CEN 51 


Early B 


14187 


9.79 


3.5 












<29.0 


<29.5 


CEN 24 


B2 


15126 


8.42 


3.5 












<29.0 


<29.5 


CEN 49 


B2 


14649 


9.72 


3.5 












<29.0 


<29.5 


CEN 48 


B2 


13124 


10.98 


3.5 












<29.0 


<29.5 


CEN 90 


B2 


14345 


11.51 


3.5 


406 


0.1 


79 e 


22.2 


1.8 


30.79 


31.39 


CEN 85 


B2 


6104 


10.62 


3.5 


466 


0.0 


538 


22.3 


4.8 


31.86 


32.19 


CEN 84 


B2 


6109 


11.11 


3.5 


495 


0.2 


8 


21: 


2: 


28.8: 


29.1: 


CEN 96 


B2 


14492 


10.81 


3.5 












<29.0 


<29.5 


CEN 89 


B2 


7390 


11.33 


3.5 












<29.0 


<29.5 


CEN 99 


B2 


9160 


11.58 


3.5 


638 


0.3 


88 


21.8 


2.2 


30.67 


31.12 


CEN 17 


B3 


12345 


12.33 


3.2 


182 


0.2 


12 e 


21.9 


0.7 


29.06 


30.47 


CEN 93 


B3 


14349 


9.59 


3.2 












<29.0 


<29.5 


CEN 35 


B3 


14179 


10.42 


3.2 


186 


0.1 


30 


22.1 


3.0 


30.42 


30.81 


CEN 75 


B3 


23809 


11.99 


3.2 












<29.0 


<29.5 


CEN 94 


B3 


12342 


11.47 


3.2 












<29.0 


<29.5 


CEN 91 


B3 


14141 


11.25 


3.2 


312 


0.1 


86 


22.1 


3.7 


30.94 


31.29 


CEN 83 


B3 


12337 


11.28 


3.2 


350 


0.0 


309 f 


22.1 


5.0 




31.85 


CEN 14 


B3 


14893 


12.16 


3.2 












<29.0 


<29.5 


CEN 26 


B3 


13500 


10.53 


3.2 


376 


0.1 


78 


22.3 


1.9 


30.81 


31.38 


CEN 74 


B3 


27036 


11.77 


3.2 


399 


0.2 


6 


21.: 


2: 


28.7: 


29.0: 


CEN 46 


B3 


7804 


9.47 


3.2 


450 


0.8 


6 


21.: 


2: 


28.7: 


29.0: 


CEN 47 


B3 


9089 


10.80 


3.2 












<29.0 


<29.5 


CEN 65 


B3 


3592 


11.50 


3.2 












<29.0 


<29.5 


CEN 27 


B3 


7284 


9.90 


3.2 












<29.0 


<29.5 


CEN 44 


B3 


1183 


13.86 


3.2 












<29.0 


<29.5 



Note. — Column 1: This list is obtained from Appendix A of TFM03 (with IRS 15 added) whic h gives optical cross-identifications , positi ons and spectral 
types. The stars arc listed first in order of decreasing mass, and then by right ascensi on. CEN = [ Chilli, Elsaesser. & NeckeH |l98CH , OI ^ lOgura Sz Ishidal 
<1976| ). The CEN 1 04+04 binary ( J5~^| is also known as Klcinmann's Anonymous Star <KlemmamrJ|1973| ) and CEN 3 is BP -16° 4818. The B3 stars CEN 50, 
CEN 58, CEN 81, and CEN 82 are omitted because accurate positions arc not availabl e. The stars CEN 33 «• B 324, C EN 34 4> B 358, CEN 95 O B 266, 
and CEN 102 & B 305 arc omitted because they are likely to be background red giants I Hanson, Howarth. & Conti 1997); none of these sources was detected 
in X-rays. 

Column 2: Spectral types arc from Hanson, Howarth, Li Conti il997), Chini, Elsaesser. & Ncckel il980), and SIMBAD. 
Columns 3-4: Source numbers and RT-band magnitudes are from | J i anl^t^dTHoO^ 

Column 5: Bolometric l umino sities are estimated from calibrations of L boi with spectral type, I Martians et alj ^2005f) for 03-09.5 stars and 
dc .lager & Nieuwenhuiizen 11987) for B stars. No use is made of available photometry. 
Column 6: C/iandra source number, from Table^ 

Column 7: Offset (in arcscconds) between the Chandra and SIRIUS sources. Both fields arc aligned to the 2MASS astromctric frame. 

Columns 8—12: X-ray properties from Table extracted counts after background subtraction; column density and plasma energy from fits to the ACIS 
spectra (: in log N H values arc approximated from median energy i Feigclson ct al. 2005); : in kT are assumed): observed hard band luminosity (2 — 8 keV); 
inferred total band luminosity corrected for absorption (0.5 — 8 kcV). Upper limits to luminosities are matched to the fainter sources in Table f3l 

a Wc adopt the nomenclature CEN la and CEN lb for the 04+04 binary components of Klcinmann's Anonymous Star ( f5~2l which arc resolved by 
Chandra. These sources do not appear in the SIRIUS catalog due to crowding and/or saturation. The K-bsmd magnitudes shown arc estimated by splitting 
the K = 5.84 measurement from fChini fc Warga^ (199 8| ) into two equal halves. 

b This star is saturated in the SIRIUS data and docs not appear in the catalog. Wo use here position and photometry from the 2MASS catalog: 2MASS 
18203586-1615431 with ,1 = 7.97, H = 7.21 and K = 6.82. 

c This star is outside the field of view of the SIRIUS data. We use here position and photometry from the 2MASS catalog: 2MASS 18200299-1602068 with 
J=8.68, H = 8.09, and K=7.80. 

d This star is outside the field of view of the SIRIUS data. The available position is too inaccurate to reliably match to 2MASS. 
c Thc X-ray emission may be variable. 



f Thc X-ray emission is definitely variable; sec Figurc[3]for light 
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5.1. OB Star Luminosities 

The most fundamental question concerning OB X-ray emission for single stars is whether the X-ray 
production scales with the wind power, which in turn is expected to scale with the bolometric luminosity. 
Figure IT3la compares X-ray and bolometric luminosities for the cataloged M17 OB population. The X-ray 




I I i i I I I I I I I ii I - ^ • I I I I I I L 

31.0 30.S 18:20:30.0 29.5 23.0 

Right Ascension (J 2000) 



Fig. 11. — ACIS image of the central region of the NGC 6618 cluster with source sequence numbers and 
extraction regions indicated. This image is shown at high-resolution with 0.25" x 0.25" pixels. Sources #536 
and #543 (bold-italic) are the two 04 components of Kleinmann's Anonymous Star ( ^5 . 2|) . Sources #495, 
567, and 574 (bold- italic) are also known OB stars (Table [9]). Sources with underlined labels are listed in 
Table [5] as candidate new intermediate- or high-mass stars. 
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Fig. 12. — Thermal plasma models for the X-ray spectra of selected OB stars in M17. The absorbing 
column density Nh is shown in units of 10 22 cm~ 2 ; plasma temperature is shown in units of keV. Elemental 
abundances were Z = 0.3Zq unless otherwise specified. Panels a, b, and e are two-temperature models; 
model components are shown separately as dotted curves. (a,b) Spectral models for the binary components 
of Klcinmann's Anonymous Star (#536 and #543, £)5.2p . In these two-temperature models, both 04 stars 
show a typical soft plasma; both also show a muc h harder plasma com ponent most likely caused by colliding 
winds from an unresolved binary component (e.g. Pollock et al. 2005h . (c) A soft spectrum from the 06 star 
01 345 (#433). (d) A B2 star (#466) with a hard spectrum; this star lies at the center of NGC 6618 (E_Q). 
(e) A hard spectrum from the 08 star 01 352 (#732). (f) A hard spectrum from the unidentified luminous 
source #396 (gOJ). 
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luminosity shown, Lt iC , is over the full (0.5-8 keV) band, corrected for absorption. We adopt L^ i values 
here based on calibrations between spectral types and bolometric luminosities for main sequence OB stars 
( Martins et al. 20051 : de Jager fc Nieuwenhuiizerilll987t ). No attempt is made to use available photometry 
with absorption and bolometric corrections as the available measurements and reddening estimates are 
incomplete and unreliable. 



We used the ASURV survival analysis package (jlsobe et al. 1986]^] to perform statistical tests for corre- 
lation between between L t c and L^oii with proper statistical treatment (generalized Kendall's tau correlation 
test) of the available upper limits. The correlation among M17 sources was very signific ant (null hypoth- 
esis probability P < . 01%) , supporting the long-standing L x oc 10~ 7 i& ; relationship (jChlebowski et al 



1989c IBerghoefer et al.lll997l ). Most of the lower luminosity members are undetected; this group would be 
underrepresented in a plot that did not include X-ray upper limits based on a sample defined at other 
wavebands. 



- 1 Available from the Center for Astrostatistics at http://astrostatistics.psu.edu/statcodes 
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Fig. 13. — Relationship between X-ray and bolometric luminosities for M17 OB stars in Tabled Arrow 
tails represent estimated upper/lower limits; overlapping points are artificially displaced horizontally for 
readability. Boxed points represent estimated values. The dotted line represents the model L tyC = 1 x 
10~ 7 Lboi- (a) Broad-band (0.5-8 keV) X-ray luminosity corrected for absorption vs. bolometric luminosity, 
(b) M17 OB star s (sam e symbols as in panel a) shown with Pismis 24 OB stars (red pluses and red arrow) 
from IWang et all J2006h and Orion Nebula Cluster OB stars (green x 's and green arrow) from the COUP 
study |Stelzer et al-lbooih . 
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For comparison, Figure [T3b shows the M17 OB stars in context wit h the OB stars in the Pismis 24 
cluster in NGC 6357 (jWang et al.lbood ) and in the Orion Nebula Cluster (jstelzer et al.ll2005r ), with L boi for 
these sources estimated from spectral types via the same method used for the M17 sources. The red arrow 
represents a lower limit for L^d for the Wolf-Rayet + 07 binary WR93 in NGC 6357 and the green arrow 
represents an upper limit to the X-ray emission from a COU P B star. The L t . c v alue for the Orion source 
1 Ori C was replaced by a more accurate value reported by iGagne et al.l (|2005r ). This plot demonstrates 
that the M17 OB stars display a similar correlation between L tjC and L(, ;, with similar scatter, to that seen 
in the other two populations. As in the COUP sample, more scatter is seen in the M17 B stars than in the 
O stars; these B stars have emission comparable to lower mass pre-main sequence stars and sometimes that 
emission may actually come from unresolved low-mass companions rather than from the B stars themselves. 
Only one of these M17 sources is variable: the B3 star =#=350 CEN 82) shows factor of —2 variations 
during the 12 hour exposure (Figure [3]). The behavior does not resemble the common fast- rise-slow-decay 
shape seen in low- ma ss stellar flares, so it is possible that it arises from large-scale shocks in a magnetically 



channeled wind fe 



NGC 6231 bv lSana et al.l (|2006h . 



Gagne et al.l 120051 ) . Results shown in Figure [T3] are consistent with those found for 



5.2. The central 04+04 binary 



Chandra resolves fFigurc fTTj) the two luminous components of Klcinmann's Anonymous Star (jKleinmann 



19731 ). the binary 04+04 system with separation 1.8" (-3000 AU projected) and total mass M ~ 140 M at 



the top of the stellar IMF of NGC 6618. The spectra l types for these stars were originally determined by 



UBVRI photometry (jChini. Elsaesser. fe NeckelH l980) using images where the two s tars were not resolvec 



thus the identical spectral types is an assumption. lHanson. Howarth. fe Contil (|19970 obtained a NIR spec- 



trum for the unresolved pair and deduced an IR spectral type of 03-04, again with no information on the 
individual components. Given that these early O stars provide substantial ionization to power the M17 
complex and lie near the center of the NGC 6618 cluster, independent spectral types for the two components 
via improved NIR and visual spectroscopy are clearly needed. 

ACIS sources #536 and #543, the X-ray counterparts to the 04 stars, are the brightest and most 
luminous X-ray stars in the cluster with intrinsic total band (0.5-8 keV) luminosities logL tjC ~ 33.2 and 
33.3 ergs s _1 , respectively. We see no enhanced X-ray emission between the two 04 stars, implying that 
there is not a strong shock between their winds; this result is expected given their large separation. The 
light curves show that the X-ray emission from each source is constant within —10% during the 12 hour 
exposure. 

Two-temperature thermal plasma models were needed to account for each 04 star's broad X-ray spec- 
trum, and each star shows a strong hard component with emission extending beyond the 8 keV reflectivity 
limit of the Chandra mirrors (Figure [T2"la.fc). We recorded twice as many ACIS counts for source #543 as 
for #5360 In both sources, the hard component contributes substantially to the intrinsic flux. For #543, 
96% of the hard-band intrinsic flux comes from the kT = 10 keV plasma; this component contributes 26% 
of the soft-band and 52% of the full-band flux (logLt ;C — 33.0 ergs s -1 ). For #536, 71% of the hard-band 



11 Source #543 may suffer slight photon pile-up with ~ 0.3 counts per frame, but a careful analysis of spectra from annular 
regions in the PSF wings of this source did not yield substantially different spectral fit parameters; although the core might be 
slightly piled up, there are not enough counts in the wings to yield a measurable change in the fit parameters when core events 
are excluded. 
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flux is due to the kT — 13 keV component, while just 10% of the soft-band and 25% of the full-band flux 
(logL t . c ~ 32.5 ergs s _1 ) comes from this hard plasma. 

Skinner et al. ( 2002h provide a particularly cogent and helpful review of the possible causes of hard X- 



ray emission in massive stars. They describe an XMM-Newton observation of the Wolf-Rayet star WR 110, 
which displays the expected soft plasma (kT = 0.55 keV) seen in single massive stars and an additional hard 
component with a poorly constrained temperature of kT > 3 keV. They assess several models, including 
magnetically-confined wind shocks, and finally conclude that wind shocks with an unseen binary companion 
are the most likely cause of the hard component. 

Another model mentioned by Skinner et al. ( 20021 ) involves non-thermal hard X-ray emission caused 
by inverse Comp ton scattering of UV photons from the star off of relativistic charged particles created in 
the wind shocks (jChen fc Whitdll99lh . Although the WR 110 data are too poor to constrain this model, 
we can apply it to the M17 04 stars. Fitting the fainter 04 star (#536) with a power law instead of 
a hard thermal component (keeping Njj and the soft plasma component the same as in the original fit) 
yields an acceptable fit, only slightly worse than the two-temperature thermal plasma fit, with a power 
law photon index of T = 1.3. Although thi s power law slo p e is n ot well-constrained by these data, it is 
consistent with the slope of 1.5 predicted by IChen fc White! (|199ll ). The full-band intrinsic luminosity for 
th is power law compone nt is log L t . c ~ 32.4 ergs s -1 , thr ee orders of magnitude brighter than the prediction 
of IChen fc Whitel (jl99ll) . However. Iskinner et alJ (|2002h note that this luminosity is strongly dependent on 
the magnetic field strength of the star and its wind temperature and mass loss rate. Thus, inverse Compton 
scattering could be an alternative explanation to colliding winds for source #536, but extreme physical 
parameters would be required to generate such a high luminosity. 

In contrast, fitting the brighter 04 star (#543) with this soft plasma -I- power law model yields an 
unacceptable fit; this source has an Fe line that the power law model cannot reproduce (see Figure fT2b ) . 
This Fe line is strong evidence that the 10 keV plasma in source #543 has a thermal origin, most likely 
caused by colliding winds between the 04 star and an as yet undiscovered binary companion. Since the 
Fe line in source #543 argues that its hard component is thermal, it is reasonable to assume that the hard 
component seen in the other 04 star (source #536) is also thermal; it is probably just too faint for its Fe 
line to be seen in this observation. 

This implies that each 04 component of Kleinmann's Anonymous Star may itself be a massive binary 
and that the system actually consists of at least 4 massive stars. It is remarkable that these two X-ray 
sources show such similar and yet highly unusual spectra; thermal plasmas this hot are not seen even 
in Wolf-Rayet colliding-wind binaries o r in r) Carinae. Only the "7 Cas analogs" ( Smith fc Balona 20061 : 



Rakowski et al 



2006 



Smith et al.ll2004l ) show such hard X-ray spectra, yet M17's 04 stars do not show the 



variable lightcurves typical of these objects. Some complex process, likely involving both colliding winds in 
close binaries and unusual magnetic field properties, may be necessary to explain M17's 04 stars. Again, 
visual and NIR spectroscopy of these two stars is highly warranted. 

The NGC 6618 04 stars and later OB stars (Figure fT2]) thus give new examples of the extremely hot 
plasmas that can be generated in powerful O star winds. The luminous, hard emission of these stars has 
implications for the study of Galactic massive star formation in the X-ray band. Emitting ~10 315 ergs s _1 
above 5 keV, such stars can be detected in reasonable Chandra exposures anywhere in our half of the Galactic 
disk, even if they lie behind several spiral arms and are embedded in dense molecular clouds; they can even 
be studied in the Magellanic Clouds. Recently studied exa mples have been found in the MSFRs NGC 6231 
NGC 3576, Wd 1, W3, W51A, and R136 in 30 Doradus |Sana et alJbood : lTownslevll2006ci: Iskinner et al 
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20061: llbwnslev et alj|2005l l2006bh 



5.3. Possible new OB stars 



Chandra data have the potential for identifying new OB stars in MSFRs, either in the rich clusters 
illuminating Hn regions or deeply embedded in the surrounding molecular cloud. Recent Chandra studies of 
RCW 38 and NG C 6357 have identified 31 and 24 candidate new OB stars, respectively (jWolk et al.ll2006l 
Wang et alJl2006h . Such stars are bright in the NIR and already appear in IR catalogs, but cannot be easily 



distinguished from reddened field stars (e.g. red giants, which are generally X-ray faint) without X-ray study. 



Jiang et al.l (|2002f ) find about 1000 OB candidates from IR photometry, but many of those are likely 
foreground or background stars. X-ray emission should help to discriminate interlopers from M17 OB stars. 
In Table [51 we narrow these candidates to 143 X-ray selected sources that have masses > 2 M Q based on the 
NIR color-magnitude diagram in Figure [Sj Twenty-nine of these are known OB stars from Tables [5] or [5J 
but the remaining 114 sources are candidate massive stars. Two of these sources, #51 and #163 (marked 
in Figure [6]), appear to be previously unremarked highly obscured O stars with Ay ~ 40. Spectroscopy of 
these OB candidates is warranted to confirm their spectral types. 



Two alternative indicators of candidate OB stars have been employed (e.g. IWang et al.l 120061 ): high 
X-ray luminosity, and X-ray detection coupled with bright if-band magnitude. There are just four X-ray 
sources (described below) in our sample with logLt jC > 32.0 ergs s _1 that are not known O stars (Tables [8] 
or [9]) ; none of them appears in our list of candidate protostars (Table [7} . 



ACIS #187 B335) This source, located in the South Bar, appears in Table [H where it is estimated 
to be an early B star (M ~ 11 M Q ). No X-ray variability is seen. Its X-ray luminosity and spectral 
characteristics are similar to the known B2 star CEN 85 (<=> ACIS #466), so it is plausibly a new 
high-mass M17 star. 

ACIS #230 This source is located near the top of the ACIS field, west of M17 North. It shows no 
X-ray variability and its IR colors do not imply that it is massive. Its X-ray spectral parameters are 
essentially identical to #187 and #466, however, suggesting that it too may be massive. 

ACIS #396 Ranking by detected counts, this source is the second brightest X-ray source in the field, 
located in the South Bar near the Hn region interface. It is not listed as a bright NIR source, but 
IR photometry can be affected by nebular emission. Its X-ray lightcurve is classified as "definitely 
variable" but this variability is not flare-like. This hard X-ray source is equally well fit by a very 
hard thermal plasma (kT > 10 keV, Figure \12}f) or a power law with photon index T = 1.55; both 
fits give the same absorption and flux. Although its NIR colors do not indicate a massive star, this 
source's X-ray characteristics are similar to the most massive stars in M17. Its flat X-ray spectrum 
could indicate that it is a bright AGN, but its location near the ionization front argues against that 
interpretation. 

ACIS #655 (<*=> B148) This source, located near the center of NGC 6618 (unfortunately in a chip gap), 
is quite faint for most of the observation, but exhibits the strongest flare in the whole sample (see 



Grosso et al. 


2004: 


Wang et al. 


2006; 



source #655 is estimated to be an intermediate-mass star (M ~ 2 M Q ), suggesting that the flare may 
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be from an unresolved lower-mass companion. It is unlikely that source #655 is a new high-mass 
member of M17. 

The final indicator of OB candidacy that we consider — X-ray detection coupled with bright (K < 10) 
if-band magnitude — produces 28 candidates. Seventeen of these are known massive stars and two are likely 
foreground stars: #281 <^ B293 <=> CEN 7 (spectral type F8) a nd #854 O TYC 6265-1977-1 (sp ectral 



type G with large parallax). ACIS #600 (<£=> B 163) is listed by lHanson. Howarth. &: Contil (| 19971 ) as a 
candidate young stellar object. The remaining sources are ACIS #1, 9, 35, 402, 650, 698, 852, and 875, 
widely distributed across the ACIS field of view. These are all listed in Table [6] 



5.4. Protoplanetary Disks Around Higher-Mass M17 Stars 



As discussed in <JT] and IFeigelson et al.l ((2006) , X-ray surveys of young stellar populations give samples 
that are largely unbiased with respect to the presence or absence of protoplanetary disks. They thus provide 
an unusual opportunity to study disk evolution in coeval populations. Here we address this issue in a 
restricted fashion, considering only inner disks revealed through a color excess in JHK photometry ( £|2.4.2|) 
arising from their heated inner edges; such disks are present mainly in actively accreting Class I and II 
young stellar systems. Our analysis is insensitive to centrally-cleared, outer disks emitting only at longer 
wavelengths. 

The stellar population we consider here is our X-ray-selected sample of intermediate- and high-mass 
stars (Table [6] in ^2.4. 2[) . This X-ray sample is not complete even for the 1 Myr old NGC 6618 population. 
Table |6] is missing disk/envelope evolutionary classifications for some known massive stars (e.g. the 04+04 
binary) due to imperfections in the JHK photometry (i j2.4.2p . Half of the cataloged early B stars in the 
NGC 6618 cluster are undetected ($5|) in X-rays, and a larger fraction of late-B stars will be missing. Our 
shallow X-ray data miss much of the massive pop ulation just now formi ng a nd thought to possess dis ks, e.g. 
many of the high-mass C lass I sources studied by Nielbock et al. ( 2001 ) and lHoffmeister et al. ( 20061 ). M17- 
SOl (jChini et al.ll2004bl ). and CEN 92 (jChini et al.ll2005f ). Despite this incompleteness, our X-ray sample of 
higher mass stars in the coeval NGC 6618 population should be largely unbiased with respect to the presence 
of inner disks. 

Our sample thus consists of the 138 stars in Table [6] with disk/envelope evolutionary classifications. A 
remarkably low fraction (12%, or 16/138) are designated Class I or Class II, indicating AT-band excesses 
consistent with inner disks. A similar result emerged from the recent Chandra study of the rich Pismis 24 
cluster illuminating the NGC 6 357 Hii region, w here the X-ray selected NIR Class III:II ratio was ~25:1 
among intermediate-mass stars ([Wang et al.ll2006l ). Considering that the NGC 6618 cluster is thought to be 
only ~1 Myr old (2]), we infer th at the inner disks around its more massive star s evolv ed very rapidly. A 
similar conclusion was reached by Hillenbrand et al. ( 1993h and Hernandez et al. ( 2005^ in their studies of 
young stellar clusters with many disk-free intermediate- mass stars. It also sugges ts that studies of nearby 
accreting Her big Ae/Be disks with strong NIR emission ( Waters fc Waelkensl 19981 ) may characterize only a 
small fraction of the underlying population of coeval stars. Disk-free stars are often missing from IR-selected 
catalogs of pre-main sequence populations. 
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6. Summary 

Extended X-ray emission had been detected from the M 17 Hn region by the Einstein, ROSAT, and 
ASCA X-ray observatories since the 1980s, but the finding had not been reported until the Chandra X-ray 
Observatory revealed a wealth of morphological detail. Our earlier study focused on the outflowing 10 MK 
plasma attributed to shocked O star winds (TFM03), while the present study concentrates on the stellar 
population. In a 40 ks ACIS-I exposure obtained in 2002, we find 886 unresolved X-ray sources with a 
wide range of intrinsic X-ray luminosities (Lh.c — 10 28 to almost 10 33 erg s _1 in the hard 2 — 8 keV band) 
and absorptions (column densities Njj < 10 20 to 10 24 cm -2 ). Nearly 90% of the sources have counterparts 
in NIR (2MASS, SIRIUS) or MIR (Spitzer GLIMPSE) images; the remaining sources are divided between 
background contaminants and newly discovered members of the cloud population. With Chandra, the 
identification of X-ray sources with individual stars is unambiguous (median offset 0.24" from NIR sources) 
except for components of multiple systems. 

The spatial distribution of X-ray stars suggests a more complex morphology of star formation than 
is generally discussed from optical- and IR-derived surveys. In addition to the central NGC 6618 cluster 
and three well-studied embedded groupings (M17-UC1, KWO, and M17-North), we find evidence for a 
new embedded cluster (designated M17-X) located 4' (2 pc) north-northwest of NGC 6618 (Figure US*)- A 
5'-long arc of several dozen X-ray stars along the interface between the M17 Hn region and the M17-SW 
molecular cloud core is c learly traced , delin eating recently triggered star formation and confirming the similar 



structure seen in NIR bv l Jiang et al.l(|2002h . Substructures on scales of 0.1 pc are seen both within the central 
NGC 6618 cluster and in the embedded populations. These suggest that the populations are dynami c al and 
young, and that equilibrium has not been reached. Our data support the argument of I Jiang et al. (|2002h 



and others for widely distributed star formation that is not concentrated into clusters. This is seen in both 
the heavily and lightly obscured X-ray populations. 

We use our detailed knowledge of the Orion Nebula Cluster population and XLF to calibrate the observed 
XLF of the ACIS M17 field. The inferred total population of heavily obscured stars is 5000 - 7000 down to 
0.1 M Q , outnumbering the lightly obscured (Ay < 10 mag) population of 2500 — 3500 stars which are mostly 
concentrated in t he central NGC 66 18 cluster. These values are consistent with the ~3600 Class I-III young 
stars counted by Jiang et al. (2002) in their NIR survey sensitive to K < 19 mag. 



Forty percent of the ACIS sources are heavily obscured, defined here by median energy E me< n an > 
2.5 keV which is roughly equivalent to Ay > 10 mag. Concentrations of 10 — 20 sources are seen around 
each of the three well-studied star formation subregions in the field: IRS 5/UC1, the KWO, and M17-North. 
Some high-L/^, high-mass embedded sources are detected in X-rays (e.g. IRS 5S and Anon 1 near M17-UC1, 
IRS 2, CEN 31, IRS 15), but others arc not (the KWO, M17-S01, CEN 92). The situation for the most 
famous embedded objects is complicated: IRS 5N which ionized M17-UC1 is only tentatively detected with 
<5 photons; the KWO Herbig Be star is undetected but a nearby member of the KWO cluster is seen. 

A much greater number of heavily obscured stars are seen outside of the three known concentrations. 
Many are located in the North Bar and South Bar where the Hn region ionization front is propagating into 
the cloud, but others appear in more distant molecular cloud cores or distributed more-or-less uniformly 
across the field. At our current completeness limit for heavily obscured stars (roughly estimated in <JT2] a $ 
logL/i c > 30.9 ergs s _1 ), only the tip of the embedded XLF is detected. An upcoming deeper Chandra 
observation should identify hundreds more heavily obscured stars and elucidate any clustering patterns of 
the embedded stellar population. 
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All of the cataloged O stars but only half of the cataloged early B stars (B0-B3) in the field are detected 
in this short 40 ks ACIS exposure. These OB stars show intrinsic X-ray luminosities consistent with the 
long-standing L x oc 1CP 7 L^i relationship, with scatter similar to that found for the NGC 6357, NGC 6231, 
and Orion Nebula Cluster populations. 

Only the 04+04 binary components, which are resolved at 1.8" separation, are very X-ray luminous, 
emitting ^2 x 10 33 ergs s" 1 in the total band, 20 — 30% of which appears in the penetrating hard band. 
Although there is no evidence for a shock between the winds from the two resolved components, each 
component itself exhibits remarkably hard thermal plasma emission even harder than that usually seen in 
colliding wind binaries. This suggests that Kleinmann's Anonymous Star may be a system of at least 4 
massive stars, two pairs of colliding-wind binaries both with extraordinary wind and/or magnetic properties. 

Many B stars have X-ray luminosities scattered in the 29 < log if < 31 ergs s _1 range, overlapping the 
low-mass pre-main sequence XLF. Most of these stars are thus not readily distinguished from lower mass 
stars in the X-ray band and, not infrequently, the emission likely comes from a low-mass companion. Many 
likely new OB stars are found by the association of Chandra sources with NIR sources. 

Our final result concerns the frequency of inner protoplanetary disks around the X-ray selected high- 
and intermediate-mass stars. X-ray selection is effective in establishing the underlying population of these 
stars in a disk-unbiased fashion. We find that only ~ 10-15% of NGC 6618's OBA stars have disks with 
if-band excesses. 
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to Gordon Garmire, the ACIS Principal Investigator, by the National Aeronautics and Space Administration 
(NASA) through NASA Contract NAS8-38252 and Chandra Contract SV4-74018 issued by the Chandra X- 
ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf 
of NASA under Contract NAS8-03060. This publication makes use of data products from the Two Micron 
All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by NASA and the National Science Foundation. 
This work is based in part on observations made with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. This publication 
makes use of data products from the Midcourse Space Experiment. This research made use of the SIMBAD 
database and VizieR catalogue access tool, operated at CDS, Strasbourg, France. We would have been lost 
without the invaluable tools of NASA's Astrophysics Data System, and without CIAO, XSPEC, and ds9. 

Facilities: CXO (ACIS). 

A. Estimation of Matching Reliability 

Because the SIRIUS, 2MASS, and GLIMPSE catalogs have far more entries than our X-ray catalog, 
some of the associations between X-ray and IR sources reported in Table [5] are expected be false positives, 
i.e. false associations found by chance positional coincidence. To estimate the number of false positives 
expected when matching to a specific IR catalog, it is useful to view the X-ray catalog as a mixture of two 
populations: an "associated population" of X-ray sources where a "true" counterpart does exist in the IR 
catalog, and an "isolated population" of X-ray sources where no counterpart exists. 

We can easily simulate the behavior of the matching algorithm for the isolated population by applying 
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an offset between the catalogs, running the algorithm, and tabulating the results. The offset must be large 
enough to "break" all the true associations, but excessive offsets should be avoided in order to preserve the 
relative spatial distributions of the two catalogs. We interpret the matches produced by this simulation 
as "false positives." We refer to the X-ray sources which (correctly) did not produce matches as "true 
negatives." 

With slightly more effort we can simulate the behavior of the matching algorithm for the associated 
population by creating a "fake" IR counterpart for each X-ray source, running the algorithm, and tabulating 
the results. Each fake IR position is randomly offset from the corresponding X-ray position by sampling the 
appropriate offset distribution, taken to be Gaussian with variance equal to the sum of that source's X-ray 
position variance and an IR position variance chosen randomly from the IR sources reported in Table El We 
attempt to eliminate the actual associations from the IR catalog by removing the counterparts reported in 
Table [5j We also randomly offset the catalogs as described above because the "background" population in 
the IR catalog may be suppressed in the regions around the pruned sources (due to observational limitations) . 
We interpret the matches produced by this simulation as falling into two categories: "correct matches" to 
the fake entry we created, and "incorrect matches" to some other source which lies closer by chance. We 
refer to the X-ray sources which (incorrectly) did not produce matches as "false negatives." 

By repeating the simulations many times, we estimate the probabilities of each of the five match types 
described above. These probabilities, reported in Table [lOl can be used to characterize Table [5] only if we 
can estimate the true mixture ratio between the isolated and associated populations, i.e. the number of X- 
ray sources having no counterparts (thus susceptible to falsely matching the background IR sources) vs. the 
number of X-ray sources having true counterparts (thus susceptible to only those background IR sources that 
lie closer than the true match). We estimate the "associated fraction" for each catalog (/a !MASS , fAsimusi 
Saqlimpsb) ^ considering the number of X-ray sources for which no match was reported in Table [101 (410 for 
2MASS, 159 for SIRIUS, 662 for GLIMPSE). Those "negative" matches are a mixture of the true negative 
and false negative match types in the simulations (columns 4 and 5 in Table [T0|) : 
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Table 10. Simulations of matching algorithm 






Associated Population 






Isolated Population 


Catalog 


Correct Matches 


Incorrect Matches 


False Neg 




True Neg. 


False Pos. 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


2MASS 


874.2 


4.2 


7.6 




865.6 


20.4 


SIRIUS 


849.9 


28.5 


7.6 




739.2 


146.8 


GLIMPSE 


872.7 


4.6 


8.8 




861.6 


24.4 



Note. — Columns 2—4: Expected number of population (1) match types. Columns 5—6: Expected 
number of population (2) match types. 
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After solving for these associated fractions (/admass = 0.53, Jasirius = 0.79, Iaglimpse — 0.23), we 
can characterize the expected flaws in Table [5] For example, the expected number of reported matches which 
are wrong is estimated as a weighted sum of "incorrect matches" from the associated population plus "false 
positives" from the isolated population: 

4-2/a 2MASS + 20.4(1-/ A3MASS ) = 12 (2.5%) 
28.5 fAsmxvs + 146.8 (l-fA slRIUS ) = 53 (7%) 
+ 2AA(l-f AGL1MPSE ) = 20 (9%) 

Note that these are averages across the catalog; the reliability of individual matches is complex since 
it depends on both the position uncertainties (e.g. larger uncertainties lead to larger matching "footprints" 
which have larger susceptibility to spurious matches) and on the local density of IR sources (e.g. crowded 
regions are more likely to generate spurious matches). 
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